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NOTATION 


Linear  Dimensions  &  Areas 

G  Chord  of  blade 

Cw  Chord  of  wing 

Q  Diameter  of  rotor 

CL  Rotor  blade  flapping  hinge  offset 


R 
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p^r*  Distances  defined  in  Figure  3 
P  Radial  distance  along  rotor  blade 
Rotor  radius 


Wing  Area 

Non-dimensionalized  distances  in 
directions  (See  list  of  subscripts  belov) 


Xcp  Distance  of  center  of  pressure  of  wing  behind  rotor 

^  Distances  in^^~N^^  directions  (See  list  of  subscripts 

below) 

Angular  Measurement 


cX  ' 
CK^ 

(3 

6 

0 

% 


Rotor  angle  of  attack 
Wing  angle  of  attack 
Fuselage  angle  of  attack 
Angle  defined  in  Figure  8 
Angle  of  downwash  due  to  wing 
Effective  blade  pitch  at  0.75  R 
Aircraft  pitch  attitude 

Wing  incidence  measured  from  fuselage  reference  line 
Angle  between  lift  vector  with  and  without  slipstream 
(See  Figure  l) 

Blade  azimuth  position 
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FORCES  AND  MOMENTS 


H-force  coefficient 

Thrust  coefficient 

Centrifugal  force  per  blade 

Blade  sectional  profile  drag  coefficient 

Force  (See  list  of  subscripts  below) 

Lift  (See  list  of  subscripts  below) 
Moment  (See  list  of  subscripts  below) 


Ihrust 


Non-dimensiorial  velocity 


V_ 

ISTr 


COEFFICIENTS  AND  CONSTANTS 

Rotor  blade  lift  curve  slope 

Wing  lift  curve  slope  wl.thout  slipstream;  Rotor  constant 
f.lapping 

Rotor  longitudinal  flapping 
Wing  aspect  ratio 

Quartic  coefficients  (Equation  19^) 

Rotor  lateral  flapping 

Non-dimensional  coefficient  (See  list  of  subscripts  below) 
Aspect-ratio-dependent  constant 
Constant  defined  by  equation  ll6 
lnfl.0%'  ratio  defined  by  evpiation  11. 

'P'A 

Constant  defined  on  page  8 
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fcr 

V 

Vr 

w 

XL 

/A> 


Inflow  tased  on  induced  velocity  in  slipstream 
Velocity  ratio  defined  by  equation  12 
Defined  by  equation  4 
B.lade  solidity  factor  s 

TTR 

Locke '  s  Number  a 

Xi 

VELOCITIES 

Body  axis  system  longitudinal  velocity 
Induced  velocity 

Free  stream  velocity  along  flight  path 
Resultant  velocity  at  wing 
Body  axis  system  vertical  velocity 
Rotor  rotational  velocity 
Velocity  ratio 


At 

h 

nro 

Kl 

tc 

w 

I. 


MISCELLANEOUS  PARAMETERS 

Number  of  b.leides 

Symbol  for  "a  function  of" 

Moment  of  inertia  of  aircraft  about  Y-axis 
Aircraft  mass 
Number  of  rotors 
Density  of  air 

Characteristic  time  (See  Figure  I9) 

Aircraft  weight 

Blade  moment  of  Inertia  about  flapping  hinge 
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SUBSCRIPTS 


D 

F 

w 

I 

L 

o 

R 

S 

T 

W 


Aerodynamic 

Drag 

Fuselage 

Gravity 

Hub 

Inertia 

Ll.ft 

Initial 

Rotor 

Slipstream,  static 

Tail 

Wing 


DERIVATIVE  CONVENTIONS 

It  indicates  the  derivative  vith  respect  to  \/  with 

and  both  held  constant. 

2.  ^indicates  the  derivative  with  respect  to  ^  with 

only  held  constant  and  C-r  free  to  vary.  This  is  re¬ 
lated  to  derivatives  of  the  first  type  by  a  relation  of  the  form 


1<?< 


^ ^  ^ I 

"bV  ^  Z>c^  . 
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3. 


Indicates  the  derivative  with  respect  to  N/  with  both 

dLV 

and  Ot  free  to  vary,  e.g.. 


dtv 


-4. 

■^V 


\o(  -fel 


dig 

Iv  dLV 


In  order  to  find  some  flight  condition  must  be  specified. 

dLV 

For  instance^  if  it  is  specified  that  the  force  in  the  2  -  direction 
be  held  constant ^  then 
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SLMMARY' 


An  analysis  of  the  stability  characteristics  of  a  tilt-ving  VTOL  air¬ 
craft  is  made.  Ifhe  stability  derivatives  associated  with  the  wing-rotor 
combination  are  estimated  on  the  basis  of  the  assumption  that  for  the  purpose 
of  stability  analysis,,  the  forces  acting  on  the  wing  in  the  slipstream  can 
be  assumed  to  depend  on  the  vector  siun  of  the  free-stream  velocity  and  the 
induced  velocity  of  the  rotor.  !Ebe  wing  derivatives  are  presented  in  the 
form  of  general  charts  whi.ch  can  be  used  for  the  evaluation  of  the  deriva¬ 
tives  for  any  wing-rotor  combination.  Kumerical  calculations  are  made  for 
a  typical  tilt-wing  VTOL  (the  V'ertol  76)  transition.  From  these  calcu¬ 
lations  a  cixrve  l.s  obtained  showing  the  variai-ion  of  the  stability  character¬ 
istics  with  wing  tilt  ang'le .  On  the  basis  of  these  numerical  results  and 
the  general  expressions  for  the  derivatives,  an  extensive  discussion  is 
given  concerning  the  effects  of  important  physical  parameters  on  values  of 
the  derivatives  and  of  the  effect  of  these  derivatives  on  the  stability 
characteristics  of  the  aircraft. 
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nmoDUCTioK 


It  is  probably  well  to  set  forth  at  the  offset  the  ptirpose  of  this 
investigation.  It  is  not  to  analyze  the  stability  of  a  tilt-wing  VTOL 
aircraft  but  rather  to  study  the  factors  which  effect  the  stability 
characteristics.  In  order  to  obtain  useful  information  of  this  type 
it  is  necessa.ry  to  simplify  the  governing  equations  to  such  a  point 
that  the  stability  characteristics;,  e.g.  the  periods  and  damping  ratios 
or  time  constants  of  the  various  modes  of  motion  as  a  function  of  the 
tilt  angle,  can  only  be  relied  upon  to  gl.ve  a  rather  qualitative  picture 
of  the  aircrafts  response  to  a  disturbance.  Certainly,  at  the  present 
time,  it  is  necessary  to  resort  to  experimental  means  if  one  wishes  to 
obtain  accurate  information  on  the  stability  charaeteristics  of  a  pro¬ 
posed  design. 

There  are  numerous  reasons  why  the  stability  prediction  problem 
is  more  difficult  than  the  similar  problem  for  low  speed,  conventional 
aircraft.  Probably  two  of  the  most  difficult  problems, however,  are 
due  to  the  high  vel,oci.ty  of  the  slipstream  at  low  fo:n#ard  speeds  and 
to  the  rotation  of  the  wing  during  conversion.  The  effect  of  the  former 
wl.H  be  the  main  subject  of  the  present  report  and  on,Ly  li,ttle  will  be 
said  about  the  latter,  albeit  it  could  easily  prove  to  be  of  equal 
importance , 

The  scope  of  this  report  is  qul.te  limited.  Consideration  will 
be  given  only  to  VTOL  aircmft  of  the  tilt -wing  type  with  articulated 
rotors  fixed  rigidly  to  the  wing.  The  results  obtained  should,  however, 
give  some  insight  into  the  characteristics  of  similar  types  of  VTOL 
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such  as  the  tilt-ving  wi.th  conventional  propellers  or  the  fixed-wing, 
tilting-rotor  aircraft.  Al.so,  only  the  longitudinal,  dynamics  will  he 
considered  although  the  lateral  dynamics  present  similar  problems  and 
are  eq^ually  worthy  of  conside.ration. 

This  investigation  wi.ll  be  crncerned  only  with  stability  of  VTOL 
aircraft.  The  quality  sought  by  aircraft  designers,  which  is  usually 
called  "flying  qualities"  is  not  the  same  thing  as  stability  although 
it  bears  a  definite  relation  to  it.  Thus,  for  instance  an  airplane 
can  be  unstable  in  the  mathematical  sense,  l.e.  have  a  response  to  an 
initially  small  disturbance  which  grows  without  limit  as  time  Increases, 
eind  sti.ll  have  satisfactory  f. lying  qualities.  This  occurs  in  a  con¬ 
ventional  al.rc:raft  with  a  ve.ry  long  period  but  sli^tjy  unstable  mode  or 
for  the  helicopter  which  has  a  slightly  unstable  short  period  mode  near 
the  hovering  flight  condltlono  Flying  qualities  are  also  intimately 
involved  with  the  concept  of  control,  since  an.  airc:raft  could  be  so 
stable  as  to  be  Impossible  to  control.  Although  these  other  aspects  of 
aircraft  f.lyl.ng  qualities  are  equa.lly  important,  we  wil.l  be  concerned 
here  onl.y  wl.th  xnathe-matlca.!  stabil.ity. 

Consideration  of  aircraft  stability  in  the  above  sense  ultimately 
reduces  to  the  consideration  of  the  stability  of  so.lutlons  of  a  set 
of  dl.fferentlal  equations  in  the  neighborhood  of  some  equilibrium 
point.  The  degree  of  difficulty  encountered  in  determining  whether 
or  not  a  system  of  differential  equations  possesses  an  unstable  solution 
depends  on  the  nature  o.f  the  equations  obtained.  The  types  of  systems 
which  may  be  obtained  are; 
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1.  Linear  with  constant  coefficients. 

2 .  Linear  with  time  varying  coefficients . 

3 .  Non-linear 

For  systems  of  equations  of  the  first  type  the  determination  of 
the  stability  is  easily  reduced  to  an  algebraic  problem.  Although 
in  the  present  report  only  linear  systems  with  constant  coefficients 
are  obtained^  a  more  accurate  analysis  of  the  VTOL  conversion  ma¬ 
neuver  could  easily  lead  to  equations  of  either  of  the  other  types. 
Thus,  in  the  present  discussion  the  question  asked  is:  If  the  ro¬ 
tation  of  the  wing  were  stopped  at  any  point  in  the  conversion  ma¬ 
neuver  and  the  aircraft  allowed  to  come  to  equilibrium  what  would  be 
the  stability  of  this  equilibrium  point?  This  is  in  effect  assuming 
that  the  conversion  takes  p.lace  over  a  very  long  time  Interval.  In 
the  practical  VTOL  aircraft,  however,  the  conversion  time  may  be 
relatively  short,  e.g.  less  than  half  a  minute,  and  it  would  be 
reasonable  to  ask  what  the  effect  of  this  acceleration  would  be.  If 
a  constant  acceleration  is  assumed,  the  equations  which  describe  a 
conversion  will  be  linear  with  time  varying  coefficients.  General 
solutions  of  systems  of  equations  of  this  type  can  not  in  general 
be  found  although  there  are  methods  for  obtaining  approximate  so¬ 
lutions  if  the  coefficients  do  not  vary  too  rapidly. 

In  general  it  is  desirable  to  have  more  information  than  just 
knowing  whether  the  motion  is  stable  or  not.  Although  it  is  desira¬ 
ble  to  have  explicit  expressions  for  frequency  and  damping  in  terms 
of  the  physical  parameters  of  the  aircraft,  this  is  seldom  possible 
and  even  if  possible  may  lead  to  expressions  of  such  complexity  that 
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the  dependence  on  a  particular  parameter  is  difficult  to  see.  In 
the  case  of  linear  eq.uations  with  constant  coefficients,  information 
of  this  type  is  easily  obtained  hy  means  of  the  root  locus  method. 
Using  this  method  the  movement  of  the  roots  of  the  characteristic 
equation  in  the  complex  plBuae  with  variation  of  any  desired  parameter 
is  easily  found. 

The  main  problem  in  determining  the  stability  characteristics 
for  VTOL  aircraft  is  the  prediction  of  the  forces  and  moments  acting 
on  that  section  of  the  aircraft  which  is  in  the  slipstream  of  the 
rotor  or  propeller  and  in  particular  those  acting  on  the  wing.  In 
most  VTOL  with  tilting  rotors,  nearly  the  entire  wing  is  in  the 
slipstream  as  well  as  part  of  the  fuselage  Euid  the  tail  over  a  large 
part  of  the  conversion. 

The  problem  of  predicting  the  forces  and  moments  on  aaa  airfoil 
in  a  slipstream  has  been  considered  by  a  number  of  authors  (Ref.  1, 

2,  3,  5  ). 

Koning  considered  the  effect  of  the  propeller  slipstream  on  the 
wing  forces  and  moments.  His  theory  is  however,  valid  when  the  in¬ 
duced  velocity  of  the  propeller  is  small  compared  to  the  free  stream 
velocity  and  thus  can  not  be  of  any  use  for  VTOL  alrcretft  where  the 
free  stream  velocity  may  be  zero  or  very  small  and  the  Induced  ve¬ 
locity  quite  .large. 

The  more  recent  analysis  of  Rethorst  was  particularly  carried 
out  for  VTOL  aircraft  and  is  hence  not  limited  by  velocity  ratio. 
These  calculations  were  carried  out  for  a  large  number  of  wing  plan- 
forms  (Ref.  6)  and  may  be  useful  in  se,lecting  optimum  planforms  for 
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a  specific  design.  An  adequate  comparison  of  these  results  vith  ex¬ 
periment  is  not  possible  due  to  the  lack  of  experimental  work  in  this 
area.  Some  comparison  is  made  in  Reference  21. 

The  experimental  work  of  Stuper  is  quite  extensive  but  was  un¬ 
fortunately  done  before  the  interest  in  VTOL  aircraft  8Uid  hence  does 
not  include  tests  at  high  induced  velocities. 

In  the  experimental  work  of  Brenkmann  it  is  interesting  to  note 
the  effect  of  the  slipstream  rotation.  For  the  condition  of  small 
angle  of  attack  and  large  induced  velocity  the  lift  coefficient  be¬ 
comes  negative  over  part  of  the  span.  It  is  concluded  by  Brenkmann 
that  this  effect  on  the  total  lift  of  the  span  is  small  since  the  in¬ 
crease  over  one  half  of  the  span  and  the  decrease  over  the  other  aver¬ 
age  out  to  that  given  by  an  average  angle  of  attack.  This  conclusion 
is  probably  valid  as  long  as  the  wing  is  not  stalled.  If  part  of  the 
wing  is  stalled,  the  actual  lift  may  be  considerably  less  than  that 
based  on  an  average  angle  of  attack.  Since  present  tilt-wing  VTOL 
aircraft  often  operate  with  the  wing  at  very  high  angles  of  attack, 
this  could  lead  to  considerable  error  in  computing  the  lift  on  the 
wing.  See  Reference  20  for  a  more  complete  discussion  of  this  effect. 
Brenkmann  also  found  that  the  slipstream  over  the  entire  wing  had  a 
destalllng  effect  which  gave  a  significant  additional  lift. 

In  order  to  make  a  qualitative  analysis  of  the  effects  of  the 
slipstream  on  the  stability  characteristic  ^  it  is  almost  essential  to 
have  a  fairly  simple  expression  for  the  lift  on  the  wing  in  the  slip¬ 
stream  which  accounts  for  the  most  Important  effects.  In  order  to  ob¬ 
tain  an  expression  of  sufficient  simplicity  to  be  useful  for  stability 
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analysis  purposes,  it  is  necessary  to  neglect  many  of  the  effects 
vhich  might  be  important. 

For  the  ptirpose  of  the  present  investigation  the  following  method 
of  estimating  the  wing  forces  and  moments  and  their  derivatives  will 
be  employed: 

We  assume  that  the  wing  forces  are  dependent  on  an  "effective" 
angle  of  attack  and  a  "resultant"  velocity  which  are  determined  by 
the  vector  sum  of  the  free  stream  and  rotor  induced  velocities.  The 
value  for  the  induced  velocity  to  be  used  is  not  that  at  the  rotor  disc 
but  somewhat  greater  than  this.  Simple  momentum  theory  shows  that  the 
Induced  velocity  far  downstreeim  must  be  twice  that  at  the  rotor  disc , 
and  in  fact  experimental  evidence  (Ref.  7)  shows  that  a  velocity  nearly 
twice  that  at  the  rotor  disc  is  reached  within  a  very  short  (less  than 
one  diameter)  distance  behind  it.  Thus  in  computing  the  general  ex¬ 
pressions  for  the  resultant  velocity  and  its  derivatives  we  take  the 
induced  velocity  to  be  (\j~  ,  where  is  a  constant  between 

one  and  two  and  v  is  the  Induced  velocity  at  the  rotor  disc. 

The  lift  curve  slope  to  be  used  in  computing  the  lift  is  obtained 
from  the  work  of  Smelt  and  Davies  (Ref.  8).  These  investigators  derive 
a  correction  to  the  power  off  lift  curve'  slope  which  is  dependent  on 
the  thrust  coefficient  and  the  aspect  ratio  of  that  portion  of  the  wing 
which  is  covered  by  the  slipstream.  This  correction  is  based  on  the 
followl.ng  considerations;  Since  the  effect  of  the  rotor  will  be  to 
increase  the  velocity  over  that  part  of  the  wing  in  its  wake,  the  lift 
will  be  Increased  over  this  portion  of  the  wing  and  hence  there  will 
be  a  spanwise  change  in  the  lift  at  the  outer  edge  of  the  rotor  which 
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will  give  rise  to  a  trailing  vortex.  The  effect  of  this  vortex  will 
be  to  decrease  the  angle  of  attack  of  the  airfoil.  If  the  aspect 
ratio  of  that  portion  of  the  wing  in  the  slipstream  is  very  large  the 
effect  of  the  trailing  vortices  will  be  negligible  and  the  lift  will 
simply  be  proportional  to  the  resultant  dynamic  pressure  and  the 
effective  angle  of  attack  of  the  airfoil.  For  smaller  values  of  the 
aspect  ratio,  however,  the  trailing  vortices  will  decrease  the  angle 
of  attack,  the  magnitude  of  this  reduction  depending  on  the  strength 

of  the  vortex  which  will  be  dependent  on  the  magnitude  of  the  spanwise 

0 

discontinuity  in  the  lift  which  is  in  turn  proportional  to  the  induced 
velocity  or  the  thrust  coefficient.  Smelt  amd  Davies  simply  compute 
the  two  limiting  cases  and  include  a  constant  which  depends  on  the 
aspect  ratio  and  is  to  be  found  by  some  experimental  means.  Using  the 
procedure  outlined  above,  the  lift  on  the  wing  can  then  be  expressed  as. 


L  -  5  Vf?  o-c-’X'  t  I  K  ( -i)]  (1) 

where  the  above  symbols  have  the  following  meaning: 


“  V  \T 


iVsi  !X?| 

V 

v  is  the  free  stream  velocity  and  NT  the  Induced  velocity  (See 
Fig.  1). 

O-o  -slope  of  the  lift  curve  of  the  vlng  section  without  the 


(2) 

(3) 
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I 


presence  of  the  slipstream. 

S  -  area  of  wing  in  slipstream. 


Xep _ 

yi^/4  t  ^cp 


) 


(4) 


where  X  dp  -distance  of  center  of  pressure  of  wing  behind  rotors. 

aspect -ratio -dependent  constant  mentioned  above  having  value 


between  0  (  /A  f?  c>0  )  and  1  (  =0  ) . 

X  I 

/\  as  mentioned  above  has  value  between  1  and  2. 

The  particular  values  of  these  two  constants  in  Ein  analysis  are 
to  be  determined  in  one  of  the  following  ways: 

1.  Experiment;  either  from  models  or  full  scale  aircraft  of  the 
configuration  being  analyzed  or  by  extrapolation  from  experimental 
data  obtained  from  similar  wing-rotor  arrangements. 

2.  Rough  estimates  based  on  physical  reasoning,  e.g.  if  on  the 

particular  aircraft  under  consideration,  the  wing  center  of  pressure 
was  greater  than  one  rotor  diameter  behind  the  rotor,  it  would  be 
reasonable  to  assume  that  could  be  considered  to  have  the  value 

2. 

3.  Choose  values  so  as  to  give  a  best  possible  fit  to  more 
accurate  theoretical  calculations  if  available. 

The  value  of  HJ"  ,  the  velocity  induced  at  the  rotor  disc,  is  ob¬ 
tained  by  equating  the  thrust  of  the  rotor  to  the  change  of  momentum 
of  the  slipstream,  i.e. 

T=  (3fV3-)(  pS  V(,  )  (5) 
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It  is  assumed  here  that  the  ma,ss  flow  is  that  flowing  through  a 
stream  tube  of  the  same  diameter  as  the  rotor  disc  at  the  resultant 
velocity. 

The  process  of  taking  the  required  derivatives  of  the  resultant 
velocity  based  on  the  above  assumptions  becomes  quite  laborious.  How¬ 
ever,  it  was  found  that  they  could  be  computed  in  a  general  manner 
and  expressed  in  the  form  of  charts  so  that  their  values  can  be  found 
immediately  if  all  the  trim  conditions  are  known. 

An  approximate  expression  for  the  resultant  velocity  can  be  found, 
which  under  certain  conditions  introduces  only  a  small  amount  of  error, 
by  assuming  that  the  thrust  can  be  equated  to  the  change  in  momentum 
of  the  fluid  based  on  only  the  component  of  the  freestream  velocity 
\dilch  is  normal  to  the  rotor  plane.  This  expression  leads  to  large 
errors  when  the  angle  of  attack  of  the  rotor  plane  is  small  at  appreci¬ 
able  forward  velocities.  Also  the  error  involved  increases  in  inverse 
proportion  to  the  disc  loading.  Thus  the  assumption  leads  to  the 
greatest  error  for  "helicopter-like”  configurations  and  to  the  least 
error  for  "airplane -like"  configurations. 

The  remainder  to  this  report  will  be  conducted  in  the  following 
manner: 

First,  we  will  derive  the  controls -fixed,  small  disturbance  e- 
quatlons  for  the  longitudina.1  modes  of  motion.  These  equations  will 
be  applied  to  a  samp.l.e  aircraft  (Vertol  model  j6)  and  the  roots  of 
the  characteristic  equation  obtained  for  a  number  of  tilt  angles 
through  the  conversion.  The  effect  of  the  various  stability  deriva¬ 
tives  on  these  roots  will  be  determined  by  means  of  the  root  locus 
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method.  Then  all  of  the  Important  derivatives  will  he  investigated 
to  see  the  relative  magnitudes  of  the  contributions  of  the  various 
parts  of  the  aircraft  such  as  the  tail,  rotor,  wing,  etc.  and  how  these 
are  Influenced  by  changes  in  the  various  physical  parameters  of  the 
aircraft  such  as  the  lift-curve  slope,  center  of  gravity  location, 
tail  size,  etc. 

Finally,  an  attempt  will  be  made  to  draw  some  general  conclusions 
from  the  Important  results  of  the  investigation  and  to  outline  some  of 
the  Important  areas  which  will  have  to  be  investigated  further  either 
analytically  or  experimentally  in  order  to  obtain  a  sufficient  under¬ 
standing  of  VTOL  dynamics  to  enable  the  designer  to  proceed  system¬ 
atically  to  the  optimum  design. 

The  dynamics  of  VTOL  aircraft  have  been  considered  previously  in 
References  l8,  19,  and  20.  The  approach  in  I9  and  20  is  essentially 
experimental.  The  einalysis  given  in  Reference  I8  is  somewhat  similar 
to  that  given  here  although  little  consideration  was  given  in  that 
report  to  the  effects  of  the  various  physical  parameters  on  the  dy¬ 
namic  characteristics  of  the  aircraft. 
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ANALYSIS 


Assumptions 

We  now  proceed  to  develop  equations  for  the  longitudinal  motion 
of  a  tllt-wlng  VTOL  aircraft.  As  the  general  development  of  small 
disturbance  equations  of  aircraft  has  been  given  many  times  before 
(See  e.g.  Ref.  9)^  we  will  be  mainly  concerned  here  with  those  parts 
of  the  development  which  differ  because  of  the  presence  of  the  rotors 
and  the  tilting  of  the  wing.  As  many  simplifying  assumptions  as 
possible  will  be  made  so  that  the  Important  effects  of  the  rotors  on 
the  stability  characteristics  can  be  seen  without  being  obscured  by  the 
multitude  of  small  effects  which  would  have  to  be  Included  to  obtain 
a  precise  description  of  the  aircraft's  dynamic  behavior.  Thus  we 
make  the  following  assumptions; 

1.  Airframe  Is  a  rigid  body. 

2.  There  Is  no  coupling  between  the  lateral  and  longitudinal 
modes  of  motion.  This  Is  normally  true  for  conventional  aircraft  If 
all  disturbances  are  sufficiently  small,  the  aircraft  Is  Initially 
In  level  flight,  and  any  gyroscopic  effects  due  to  the  engines  are 
small  enough  to  be  neglected.  The  only  other  forces  acting  on  the 
VTOL  which  might  give  a  coupling  are  the  lateral  rotor  forces  and 
these  cancel  for  two  rotors  which  rotate  In  opposite  directions. 

3.  All  disturbance  quantities  are  small  enough  that  their  products 
and  squares  can  be  neglected. 

4.  Aircraft  performs  a  constant  altitude  conversion  and  Initially 
all  velocity  components  are  zero  except  along  the  X-axls . 
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The  'Wins  tilt  ang,Ie  d.s  fixed  dtird.ng  the  axialysi.s  as  well  as 
all  other  contro.l.s.  This  asstm).ption  could  very  well  he  untenable  ex¬ 


cept  when  the  conversion  takes  place  over  a  long  period  of  time. 
6,  The  rotor  satisfies  the  following  relations; 
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where  AJ  is  the  Induced  velocity  at  the  rotor;  T”  i.s  the  total  thrust 
produced  by  both  rotors;  N  the  number  of  rotors;  f?  the  radius  of  the 
rotor;  ,D-,  the  angular  speed  of  rotation  of  the  rotor;  ©  the  col¬ 
lective  pitch  at  0.75  R  ;  S  the  profile  drag  coefficient  of  the 
blades  and  Oj  ,  Cm  ,  X  defined  by: 
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where  V  is  the  free  stream  velocity  and  o(  is  defined  in  Figure  1. 
The  directions  of  T  and  M  are  also  defined  there.  For  the  deri¬ 
vation  of  the  expressions  for  Oij-  and  On  see  References  10  and 
11.  A  number  of  terms  were  dropped  from  the  expressions  given  in 
these  books  due  to  the  fact  that  J.s  much  less  than  one  for  most 

tilt-wing  configurations.  It  is  also  assumed  that  is  constant 

in  all  of  the  above  expressions. 

7.  All  force  and  moment  derivatives  are  constant  at  each  trim 
point . 

8.  Effects  of  rotor  on  wing  forces  can  be  described  by  the  ex¬ 
pressions  given  in  the  first  section  of  this  report.  (Expressions  1, 

2,  and  3). 

Non-dimensionallzation 

Non-dimensionalizing  the  eq,uations  of  motion  for  this  type  of 
VTOL  aircraft  does  not  J.ead  to  the  degree  of  simplification  found  in 
the  case  of  conventiona.1  aircraft.  This  is  due  to  the  fact  that  the 
rotor  forces  are  proportional  to  the  square  of  the  rotor  tip  speed 
rather  than  the  dynamic  pressure.  When  the  force  equations  are  di¬ 
vided  through  by  the  wing  area  times  the  dynaml.c  pressure,  the  wing 
and  tall  forces  are  simplified  to  the  usual  lift  and  drag  coefficients  j 


25 


the  rotor  forces,  however,  are  preceded  hy  the  factor  ) 

which  remains  nearly  constant  through  a  transition  since  -Q-R,  is 
constant  and  the  resultant  velocity  does  not  vary  much.  All  of  the 
force  terms  are  never -the -less  made  Independent  of  the  size  of  the  air¬ 
craft  which  is  probably  the  most  important  reason  for  the  non-dimension- 
alization.  The  moment  egtiation  was  also  divided  by  the  wing  chord. 
Equations  of  Motion 

A  body  axis  system  will  be  used  in  -which  the  fuselage  reference 
line  is  initJ.ally  horizontal  and  pointing  into  the  wind.  The  varia¬ 
bles,  however,  will  be  expressed  in  terms  of  the  familiar  wind  axes 
variables,  i.e.  \/  ,  instead  of  U.  ,  UT  ,  ©4'  Where  V 

is  the  magnitude  of  the  free  stream  velocity,  is  the  angle  be¬ 

tween  the  fuselage  reference  line  and  the  free  stream,  %  is  the 
angle  between  the  horizontal  and  the  fuselage  reference  line,  UL  and 
tjL.r  are  the  velocity  components  along  the  body  axes.  For  hovering 
flight  OL  and  UCT  must  be  retained.  The  relations  between  these 
■variables  are : 

\,JL  ~ 


Oj  -  V  ^  -f 


(13) 


The  inertia  forces  for  the  longitudinal  motion  only  are  then: 
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M^j  =  -  I<^  S-f 


t 


(16) 


For  the  derivation  of  these  expressions  see  e.g.  Reference  9* 

Using  the  relation  =  90+  G\  —  tw  (  C>(  ^  Lvv  axe  defined 
in  Figure  l),  we  obtain: 

F/  j  =  ~rn  Jj-  V fLKjrs.  Ka-  tw)-  Vce^  -  twX<^ "  (17) 

^7.-Z-  t  Vcevi-  (o^-Lw)-^/>U*Y^  C«3^-JUw)(^“  ©f)]  (18) 


Hu  j  =  -  9^  L 


a  *  '  T 

o  • 

^  since  Lsa,'  is  constant  at  each  trim,  point. 

The  forces  acting  on  the  airciuft  due  to  gravity  are: 


(19) 


K.j  = 

N's/  y!L*jr> 

(20) 

11 

(21) 

M 

H ' 

=  O 

(22) 

.From  Figure  2  we  see  that  the  aerodynamic  forces  are: 

F^s  +  -V  F^-r  +-  Fxw  i" 

(23) 

Ft-go^ 

Fz^;  +F^p>+  Fzj  -+  Fz-vwi"  Fzp 

(24) 
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'“^pFxp  ■"  Z.pF  i"  +  ^^0l,^!,  +  ^^aJf::T 


(The  subscripts  S  ,  R  ,T  ^W,  F"  indicate  sllpsteam. 

rotor. 

tail; 

wing,  and  fuselage  respectively;  see 

Figure  2)  where. 
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The  distances  Xr  constant  hut  depend 

on  the  angle  of  tilt  Uw  .  They  are  given  by  the  expressions: 


X/^  = 

+-  (  r  Hr  CuCi^-  t  w 

(36) 

(r  ■V-ip4  ja^)  Cjbd^L-W 

:  (.37) 

i"  ^  \a/ 

(38) 

II 

2.c]l  (—04-  -5^  •+  ^  cl.^  /^-kX\  V^va/ 

(39) 

where  X<j 

Figure  3-  These  distances  have  also  been  computed  for  various  tilt 
angles  and  are  shown  plotted  as  a  function  of  LW  in  Figure  4. 

Thus  the  sum  of  the  forces  said  moments  acting  on  the  aircraft 
may  be  expressed: 

Fxxi-  Px  f}  i"  V"x  T- 1  Fxw 

^ZT  t  Fz  t  '■V'  Fkf?  -t  r  Fz-W  4- 

M  ~  FI  ~  ”  ^O-F y.  5  -  Xp  Fx.  R  “  Fx  R 

^X.-i|-  Fx-p  “  ^yFy-y"  Xo^Fxw”  XooFvw 

+  XpFz.p  t  ZfrFyp-4-  McF  Ho^Ct 


Zc: 


lo,  J 


are  constants  and  are  shown  in 
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Since  ve  are  considering  only  the  longitudinal  equations  of  motion 
with  controls  fixed,  we  have  three  degrees  of  freedom.  The  three  vari¬ 
ables  will  be  taken  as  V  ,  p(  ,  .  As  already  noted,  derivatives 

with  respect  to  oK  are  the  same  as  with  respect  to  Q<4:-  . 

If  we  expand  all  the  forces  in  Taylor  series,  retain  only  the 
first  order  terms,  neglect  products  of  the  perturbation  quantities, 
and  write  for  simplicity  the  variables  themselves  for  their  pertur¬ 
bations,  we  obtain  the  equations  in  the  form  given  below.  The  process 
of  taking  the  first  term  of  the  Taylor  series  for  the  inertia  and 
gravity  forces  along  with  the  aerodynamic  forces  eliminates  the  necessi¬ 
ty  of  subtracting  the  steady  state  solution  from  the  perturbation  e- 
quatlons  and  also  gives  rise  to  the  convenient  symbols  such  as  C'm0 
for  the  non-dimensional  form  of  the  moment  of  inertia  about  the  Y-axis. 
The  derivatives  are  constant  for  a  given  trim  condition  but  vary  with 
the  variable  uw  •  The  perturbation  equations  are  thus  given  by: 


^  V  ♦-  ^  V  -V  ©£  +  ©£  =  O  ikk) 

hot  hQjf  '  ^0f  ' 

Evaluation  of  stabilitv  derivatives 

Due  to  the  Influence  of  the  induced  velocity  on  the  magnitude 
and  direction  of  the  resultant  velocity  at  the  wing  section  in  the 
slipstream,  most  of  the  stability  derivatives  will  be  functions  of 
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the  rotor  thrust  T  .  In  order  to  sepeirate  the  effect  of  the  change 
due  to  rotor  thrust  mriation  from  that  due  to  variation  vith  the 
thrust  held  constant,  the  derivatives  will  be  taken  in  the  following 
manner; 
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is  the  thrust  coefficient  defined  by. 


T-  p  N  'Tr  f?^  CxiR)^  Cr 


(h8) 


and  the  subscripts  beside  the  vertical  line  indicate 


variables  held  constant  while  taking  the  derivative. 
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manner: 
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which  is  sufficiently  accurate  as  long  as  “ 

compared  to  unity;  "tdiich  is  nearly  always  true  for 
definition  of  is. 


V  ■ 

is  small 

VTOL  aircraft .  The 
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In  a  similar  manner. 
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The  derivatives  cam  be  calculated  and  plotted  in  a 

general  manner  so  that  they  may  be  applied  to  any  rotor  at  any  oper¬ 
ating  condltion-o  This  has  been  done  in  Reference  12.  Expressions 
for  emd  plots  of  these  derivatives  are  given  in  Appendix  A. 

Thus  we  see  that  expressions  for  the  derivatives  of  the 

> 

-forces  and  Y-moments  with  respect  to  the  three  variables  , 
V  ,  0-{-  are  needed.  These  will  be  discussed  in  two  sections; 

first  the  static  derivatives  and  then  the  dynamic  derivatives. 

Static  derivatives 

The  static  derivatives  which  are  reipaired  are  the  derivatives 
of  the  forces  and  Fz  and  the  moment  M  with  respect  to 

velocity^  angle  of  attack  and  pitch  angle.  These  derivatives  in 
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(  Qj^  is  .ho.ld  con=!t-«'f-.  In  pfl.ch  o:f  j-.b.^"^''-:  dt..rl'7n.r-i,Tes  alsO;,  but  thl.s 
is  not  indicated  by  a  subscript,  since  no  confusion  is  .like.l.y  to  result 
over  this  point) 
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(58) 


The  derivatives  of  the  moment  about  the  Y-axis  with  respect  to 
V/  and  C?(  are  given  below.  The  lower  case  letters  ,  Z-  , 
are  the  various  distances  given  previously  non-dimenslonallzed  with 
respect  to  the  wing  chore. . 


(60) 
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We  must  now  find  expressions  for  the  derivatives  of  the  forces 
acting  on  the  various  components  of  the  aircraft  with  respect  to  the 
free  stream  velocity  and  angle  of  attack.  We  note  that  the  deriva¬ 
tives  of  the  inertia  and  gravity  forces  are  zero  since  each  term  of 
the  inertia  forces  contains  a  perturbation  quantity  which  is  zero 
vhen  evaluated  at  the  trim  condition  and  the  gravity  forces  do  not 
contain  the  variables  <jK  and  V  .  The  fuselage  force  derivatives 
can  be  carried  no  further  since  there  is  no  theoretical  method  availa¬ 
ble  for  treating  the  odd  shaped  fuselages  which  are  likely  to  be 
present  on  VTOL  aircraft.  These  will  normally  have  to  be  found  by  a 
wind  tunnel  test  on  a  model  unless  the  fuselage  is  similar  to  that  of 
a  conventional  aircraft,  in  which  case  it  can  be  estimated  by  euialyti- 
cal  means.  However,  if  part  of  the  fuselage  is  in  the  wake  of  the 
rotor,  any  type  of  analytical  estimate  will  probably  be  very  difficult. 
The  expressions  for  the  various  forces  acting  on  the  aircraft  in  terms 
of  the  appropriate  velocities  and  coefficients  are: 


Ds  -  X  P  ^ 

(61) 

t-s  *  P  S 

(62) 

Dvw  ”  XP 

(63) 

P 

(64) 

p  C  Q-f 

(65) 

k-T  ”  f  ^  ‘-T 

(66) 

T=  ^  Tr  K"-C-nR)^  Ct 

(6?) 

H  -  ^TrR^C-£LR)'‘CH 

(68) 
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Substituting  these  expressions  into  the  equations  and  carrying 
out  the  differentiation,  we  easily  obtain  the  following  expressions 
for  the  required  derivatives: 

-  [%  ^  U] 

) 


^r1  =  _^INJ  (^f  1 

M  s  V.  vr  /  L  s V  'i*  av  j 


I  S-l  "  ^l]  (7!) 

^1--  (« 
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OV 

^)Cscw\ 
dv  l« 

I' 

bCzwi 
bd  I 

^)V  M 


bC^i 
&v  w 


^CxtI 

Iv 


*  -  §w  V  euw  (77) 

^  Vjf 

r“  ■  (79) 
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(84) 


Co^6 


The  only  point  -which  might  need  some  explanation  in  obtaining 
these  expressions  is  the  evaluation  of  the  derivatives  j 

The  angle  0(~r  is  obtained  from  Figure  2  and  is, 

C^  Y  ~  So  -*r  0(  ~~  -t*  L  f  ~  ^ 

_  ae  .  -  ae  ^ 
w  w 

aQ(T  _  I  —  a^  _  I d6  a^ ' 

2}d  ^  dcX 

since  the  downwash  at  the  -tail  is  a  function  of  the  angle  of  at-fcack 
of  the  wing.  The  angle  of  attack  of  the  -wing  is. 


Thus, 


The  only  derivatives  which  -we  have  left  to  evaluate  in  terms  of 
known  physical  parameters  are  derivatives  of  Ch  ,  Vft  ,  4^  , 

,  C.|_s  ,  respect  to  (p(  ,  \/  .  We  will  now  con¬ 

sider  the  e-valuation  of  these  derivatives. 


(85) 

(86) 
(87) 


So  o<  - 

(88) 

'ho(' _ 

(89) 

d(?('  _ 

'-Itlv 

(90) 

gives. 

_ 

av  - 

1 

do(‘  av  \<^ 

(91) 

dd 

1 

< 

(92) 
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In  order  to  find  the 


c 


c- 


Cm  derivatives  ve  need  an  expression  for 
I  yU  f  «  Itiis  is  obtained  from  the  ex- 


H  in  terms  of 

pression  given  previously  by  substituting  the  expressions  for  Oh ,  b,  , 
and  Q/q  given  in  Reference  11  into  it.  Ihe  result  is: 

A*)]  '») 

Taking  the  derivative  of  this  expression  wl.th  respect  to  , 

Ct  ,  and  neglecting  powers  of  greater  than  or  equal  to  two 

as  compared  with  unity,  we  obtain  the  following: 

“  (9l() 

h  ^  ->*'  I  +  ^ 

(j:^— — ^  ^  3.b3(p^  J 

f  '8) 

(95) 

^  ‘feu.  ■  A’*)''  [&-*%,('-  -> 
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The  derivatives 


dC^  1  and  ^Cv-t  1  ,  are  then  obtained  from. 


dCM 


BGrL 


(97) 


^Ch  I  =  ^Cv-4  \  J_  "5  Ch  I  &Gr  1 

b<^  w  ^  ^Gr'f^^.V  e)c/  Iv 


(98) 


The  lie^ivatives  v«re  given  previously. 

The  values  of  the  resultant  velocity  Vf^  and  the  angle  between 
the  resultant  eind  free  stream  velocities,  ,  depend  upon  the  value 
of  the  induced  velocity  of  the  rotor.  Simple  momentum  considerations 
show  that  the  Induced  velocity  far  downstream  approaches  a  value  twice 
that  at  the  rotor  plane  in  the  absence  of  dissipation.  Since  experi¬ 
mental  evidence  shaw  (Saa  fief.  7)  the  Induced  vaikocity  apfonoachM 

a  value  twice  that  at  the  disc  at  a  rather  short  distance  behind  the 
disc,  using  the  induced  velocity  at  the  disc  to  compute  the  resultant 
velocity  acting  on  an  airfoil  behind  the  disc  could  lead  to  consider¬ 
able  error.  In  order  to  account  for  this  effect  we  will  use  for  the 
Induced  velocity  the  value  ^5"  where  ^  may  have  any  value 
between  1  and  2  to  be  chosen  later  and  lyy  is  the  value  of  the  in¬ 
duced  velocity  at  the  disc  which  is  given  by: 


T 


TTR^  SIR 


'(99) 


ko 


Consideration  of  Figure  1  shows  that  the  resultant  velocity  is 


given  hy: 

which  defines  euid  ^  .  Thus, 

Vr=  hr 


(100) 


(101) 


The  relation  between  ^5  and  is, 


-HR 


Taking  the  derivatives  of  Vr  with  respect  to  CX  > 
and  Ct  ,  we  obtain  the  following  expressions: 


(102) 

^/SIR  , 


(103) 

(104) 


CbA.a  (105) 

Using  the  charts  of  Reference  12  (See  Appendix  a)  for  and 

its  derivatives,  we  may  plot  the  derivatives  of  Vr  in  general  charts 
for  a  given  value  of  .  This  has  been  done  for  the  two  extreme 

values  of  and  the  results  are  shown  in  Figure  5*  It  may  be 


ik. 

Cni?)i 


4l 


noted  that  in  some  cases  the  charts  for  the  two  different  values  of 

differ  considerably.  Using  these  charts  with  the  trim  conditions 
we  easily  evaluate  the  derivatives  |  ^  and  |  v/ 

from  the  expressions: 


(107) 


(108) 


bVp  1  r  bVR  ^  -V  1  I  .  1  ._L  (106) 

1  _  ^^1  1  'dCv[  (107) 

From  Figure  ^  may  be  expressed  as, 

yUJn  <|>  =  Cfcr^O(  (108) 

Taking  the  derivatives  of  this  expression  with  respect  to  0*i  > 
and  we  obtain  the  expressions: 

Cc64?  -C^kLo(  ^ 

+//<(  Xs -X)!!  (109) 

(^^y^c?<-Xs)(X5^^-  ^y<UinCA 

(110) 

Ik 


(111) 


The  derivatives  were  also  plotted, in  the  form  of  charts  and  are 
shown  in  Figure  6.  All  of  the  derivatives  are  functions  of  the  parame¬ 
ter,  ^  ^  and  ,  where  ^  is  the 

value  of  ^  for  V  —  O  .  Also  included  with  these  charts  are 
charts  giving  ^  ^  smd  the  derivatives  of  the 

induced  velocity  NT"  . 

Using  these  charts  and  the  trim  conditions,  one  can  easily  obtain 
t^e  derivatives  1^1^  ,  ^rom 

the  expressions; 


bVRl  =  ' 

1  ^Vr  I 

[  + 

^Vr  1 

SCt  I  1 

\ 

( 112  ^ 

1 

U,Cr 

~3^t\ 

C7<>V 

\X±^  J 

\ 

1 

bs/R  1 
^  1 

^Qr\ 

do^  >  V 

(113) 

r^>4> 

a4> 

lo<;V 

'o' J  UR 

(11^) 

< 

l( 

1 

V,Ct 

■^T 

I 

dc^\ 

Iv 

(n5) 

As  mentioned  previously,  the  lift  coefficient  of  the  wing  in  the 
slipstream  is  represented  by  the  expression, 

Cus-  Cuo(^’)[.'+  V<  (,T?Wv 

i+3 


(116) 


where 

0|_Q  is  the 

lift  curve 

slope  of  the  airfoil  for  the  power-off 

condition,  i.e.  for 

!Vr=  o 

in  the  above  expression. 

Thus, 

1  _ 

^Ol_o 

bd\  +  Cuo  (117) 

&o<  \v  w 

euid  since. 


we  obtain, 


(118) 


(119) 


An  estimate  of  the  magnitude  of  the  second  term  in  the  expression  for 
^  shows  that  it  is  much  smaller  than  the  first  term  and 
thus  we  take. 


dCus 

b<A 


^Ct-sl 


In  a  similar  way  we  find. 


The 


G 


derivatives  are  simply  given  by: 


SCdsI  « 


^(Ibs  bd‘ \  ^ 

bCos 
6d'  ^ 


(120) 


(121) 


(122) 

(123) 
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In  order  to  make  an  estimate  of  V-y  and  its  derivatives,  two 
separate  cases  must  be  considered.  The  first  is  for  the  tail  not  in 
the  slipstream.  This  case  will  be  valid  for  hovering  and  all  high 
tilt  angles  up  to  some  angle  which  will  have  to  be  determined  from 
the  geometry  of  a  particular  tail-rotor  arrangement.  Although  the 
rotor  disturbs  the  flow  somewhat  even  outside  of  the  slipstream,  the 
disturbance  is  quite  small  (See  Ref.  13),  and  can  probably  be  neg¬ 
lected.  For  this  case  then  we  have. 


V 


I 


O 


When  the  tall  is  in  the  slipstreeun  the  velocity  can  as  a  first 
approximation  be  taken  as  Vr  computed  for  ^  and  the  deriva¬ 

tives  are  those  given  previously  for  Mr  . 

We  will  now  consider  the  derivatives  of  the  rotor  pitching  moment. 

In  order  to  do  this  we  must  first  have  an  expression  for  the  pitching 
moment  acting  on  the  aircraft  due  to  the  rotor.  This  expression  is 
obtained  in  the  following  analysis: 

Figure  7  shows  the  forces  acting  at  the  flapping  hinge  which  may 
cause  a  pitching  moment  if  offset  hinges  are  present.  If  the  angle  ^ 
is  small,  the  resultsuit  force  on  the  blade  is  approximately  equal  to 
the  centrifugal  force  and  we  may  approximate  the  vertical  force  by, 

O.F  (124) 

where, 

054.  l|;  -  b»  <4^  (125) 
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The  moment  acting  at  the  hub  is  thus  approximately. 


■X  eF|3e 


(126) 


where  S  is  the  amount  of  hinge  offset.  The  pitching  moment  about 
the  Y-axis  is  then, 


-Mh 


(127) 


We  are  oxiLy  interested  in  the  value  obtained  by  averaging  this  ex¬ 
pression  over  the  angle  .  This  is  for  blades, 

,  r 

M>^o=  -  ^  ]  CF  e 


=  k  CF  e  0-, 

A 

The  centrifugal  force  is. 


pR 

CF=  a?  ~mr  dr  =  Xl?  H: 

Jo 


(128) 


(129) 


where  Ms  is  the  static  moment.  Thus, 

^  -0.=*  Ms  e  o_, 


(130) 


For  small 


-  > 

dj  =  c-cj<i,  a  C  J"  ®  ^ 


(131) 


A  more  complete  analysis  of  this  moment  valid  for  any  ^  and  |3 


is  given  in  Reference  10. 
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The  required  derivatives  are  given  by. 


^)V 

u 

^CmcF 

5V 

i  .  ^ChcpI 

U.Ct  ^C-r  K^' 

\ 

V 

(132) 

^Cmcf  I 
d<^  1 

!V 

^CmopI 

“5^  1 

_i_  3Crr-i  C.F  \ 
IV.Ct^  ^)Ct 

■aCx. 

(133) 

The  terms  on  the  right  are  easily  obtained  from  the  above  expression 
for  the  moment  as. 


.  5i(^Qi£)*Kcp  (136) 

•^0^  \i3^,v  ^  Vr'  '  dCj 


where,  Kcp  =  N  b  e.  Ms  (13Y) 

j>  Sc 

Dynamic  Derivatives 

Probably  the  most  important  dynamic  derivatives  are  those  which 
arise  due  to  a  pitch  rate,  or  the  derivatives.  Contributions 

to  these  derivatives  arise  from  the  tall  and  the  rotors.  There  is 
also  a  small  contribution  due  to  the  main  wing  but  this  is  usually 
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From  Figure  8  we  may  calculate  the  change  In  the  velocity  at  the 
rotor  plane  due  to  the  pitch  rate  01^  .  Using  simple  trigonometric 

relations  we  find, 


A-wn 


L  \xxj<\ 


(144) 


(145) 


Since  we  ai*e  ooly  interested  in  the  product  of  these  two  relations  we 
may  multiply  them  together  and  simplify  considerably  the  resulting  ex¬ 
pression.  ThM, 


}  I w*  3 

Using  the  double  angle  relations  this  is  readily  seen  to  reduce  to. 


- _ 

Expression  for  the  distances  in  the  above  expression  were  given 
earlier . 

The  expression  for  ^  in  terms  of  0.^  then  becomes, 

V/u.-ncx-  Nr-t- 

_ 


(146) 


(147) 


(148) 


I 


and; 


SUR 


Thus, 


jj  A>Mn0 

^Qf  ^  SUR 

Using  the  expression  given  in  Reference  11  for  Ch  t 


h  ' 

The  parameters  'X  , 


r ■+  6<xt  _  kkXe  ^  3>vQ-i  ^  jUia.^  _  o-ob>  ^  [ 

we  can  obtain  an  expression  for  <-* — • 

^ 

Hj  ,  M  change  with  pitch  rate.  Thus, 

r 

^Cm  ^  bC-^  'bcbi  .  ^Ch  b>M  hX 

b^  byw  b0f  bX  b&f 

Taking  the  derivatives  of  expression  I5I  and  neglecting  the  small 
terras,  we  obtain. 


(149) 


(150) 


('I51) 


(152) 


OSS 

-■  -Til  1  BLI  » 

C  Q“i 

'  3. 

OCT 

byM 

'  SK 

bCn  :r- 

CjUT 

Ta 

9^ 

t  A  ^  ®] 


fA©t 


a 


j 


(153) 

(154) 

(155) 


The  derivative  is  Just  the  rotor  lag, 

bGf 

_  bcxi  j  ^  6cXi  1  boC  ^  ^CLi 


_  dcLi  1  ^  Qo^\ 


bOji  '  6^;^ 
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(156) 


} 

1 


where 


^0-1  i 
^  I  CXi 

The  last  two  terms  of 


^Xl 


^09 


are  found  to  he  negligible 

&Gf 

fore.  From  Fig.  8  we  see  that. 


derivation  of  this  expression  see  Ref.  10 ) 


(For 

was  evaluated  be- 


^  Co(+  /s)  c^o( 

Combining  these  expressions  we  obtain, 

I  b  ^  X0  3 


52.  .;islR 

The  derivation  of  this  expression  was  also  carried  out  using  the  more 
accurate  expression  of  Ref.  10  and  the  additional  terms  obtained  were 
fotmd  to  be  quite  small. 

The  lift  at  the  tail  is  given  by, 

p  ^  X  ^  T 

When  the  aircraft  is  pitching  at  the  rate  ©f  ,  the  angle  of  attack 


of  the  tall  is, 


<^T  *  X  Of  ~  w  + 

VT 


(158) 


(159) 


(160) 


(161) 


Thus, 


Also, 


■akx  _  j_ 

^  Q-T 

^  w 

(162) 

'  , 

It 

(163) 

dGf  ^  1 

^0<y 
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We  find  the  contribution  due  to  the  tail  fan  in  the  following  manner. 


-pit-  ^  ^  ] 

From  Figure  9  ve  see  that, 

Vr  C<T  -»•  -Pt  6f  CW-r  -  f^-p 


)\t  ~ 


CilP?)T 


(164) 


(165) 


Then, 

2>fef 


L-r  Obd-Or 

(A^)t 


Xr  -  Vr  Ccsd-O^T  —  ~ 


0{?Vt 

Thus, 

ST-r  _ 

piT  R-^  iOJR)^ 

bC, 

The  pitching  moment  derivative 


JIt  CA-d.  cX-r 

MiF 


^0f 


is  computed  in  the 


following  manner: 

The  expression  for  the  moment  was  given  eeurlier  as. 


H,j,  =  Hs  n."*  <2-  OL, 


Thus, 


^^mc4=  _  ?)C>^c.P  ^Q-»  -  Ni  b<2-  (  -  \(o  \ 

5  0f  bcL,  b&c  psc.  Vv^/  k  yxiy 


bo-i 


yii/ 


(166) 


(167) 


(168) 


(169) 
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The  derivative  of  the  X-foroe  equation  with  respect  to  ep  is 
given  hy. 

These  terms  have  been  evaluated  previously.  The  derivative  of  the 
Z-force  equation  with  respect  to  9p  is  likewise, 


(170) 


- 

'd'eip 


A 


P5Vr 


,  \-  inn  C<^  -  i- w)  +-  M  1 

.  0  ©.C  n  ** 


Tof 


These  terms  were  also  evaluated  while  finding  the  moment  derivatives. 
The  derivative  Om  ©p 

From  equation  42  we  obtain. 


(171) 


~Z>  Om  _  ^  -  —  A 

b  ©f  p  S  Vr^  b‘6f  p  S  Vr^C 

where  is  the  moment  of  inertia  of  the  aircraft  about  the  Y-axis. 

This  is  not  constant  but  varies  with  the  wing  tilt  euigle  X-vv/  .  How¬ 
ever,  the  amoxint  of  variation  is  quite  small  euid  may  usually  be  neg¬ 
lected. 

Derivatives  with  respect  to  V 


bGx  « 

a 

b'C' 

P  S  Vc/ 

PSN/R? 

bC^k  _ 

a 

- 

-  a 

bs/ 

p  5 

^5  Vr 

bM 

b  </ 

a'5' 

—  v-/ 

XTi  (cX  — 


o 


~m  C-ftO.  c<  — A. 


o 


o 


(172) 


(173) 

(174) 

(175) 
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Derivatives  vith  respect  to  rate  of  change  of  angle  of  attack  ^ 

V  (•:?<  —  i'w)© 


psvg^^  aa(  “ 

P  S  Vr=» 

-4.  ^  03CT 

A 

b» 

P  s  VrJ 

-  P  5  v«? 
+  Si  Car 

b  cx 

■&Cm 

bM 

„  SiCiT 

P  ^  C 

-  7-T 


bC: 


'3Ct 


ba( 


These  derivatives  arise  due  to  the  fact  that  there  is  a  time  lag 
between  a  change  in  the  angle  ^  and  its  effect  on  the  downwash  angle 
at  the  tall  ^  .  This  time  lag  is  the  amount  of  time  req^uired  for 

the  flow  to  move  from  the  wing  to  the  tail  and  is  given  approximately  'by, 


At  - 


vr 


The  contribution  to  the  angle  of  attack  of  the 


tail  due  to  this  effect  is  then. 


^  ^  At. 
'^o(  bt 


o(  ^ 


euld  the  suigle  of  attack  of  the  tail  is, 

•+•  CX  —  l-vs/  — 


The  derivative  of  this  with  respect  to  is, 

~(dc^~r  _  i?.y 

The  derivatives  of  CjLx  given  by, 

-  -I_3 _  y  -t  CCri-tl 


p  5  vr^  i-  ~d3r 


(176) 


(177) 

(178) 


(179) 


(180) 


(181) 


(182) 
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Thus, 


V-r 

C^R)t 


CJbtft  a  -y 


Comhinlng  all  of  the  above  expressions,  we  obtain, 


(190) 


Q.~r 

5  be/ 


I^QLt  C  StVt  i*  CA4(7(t 


kiJO  ^  i- 


StVt 


bC 


b<3(T 


2T  C_^ 


(191) 


^CzLT  «  —  r  o-T  cjba  e 

Zi6(  ^  ^ 

~  St  V-r  (192) 


The  most  important  parameters  needed  in  order  to  evaluate  these 
expressions  are  the  downwash  angle  at  the  tail  €:  and  its  deriva¬ 
tive  with  respect  to  angle  of  attack  (A  .  These  derivatives  are 
very  difficult  to  evaluate  with  accuracy,  since  the  flow  at  the 

tall  is  effected  considerably  by  both  the  wing  and  the  rotor  and 

* 

the  effect  of  these  change  with  the  wing  tilt  angle  Lvs/  » 
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Evaluation  of  Coefficients  to  Stability  Quartlc 


The  non-dimensional  derivatives  developed  in  the  previous  section 
can  he  used  in  place  of  the  dimensional  derivatives  in  equations  4-3-45 
since  all  of  the  force  derivatives  were  non-dimenslonalized  with  re¬ 
spect  to  the  same  force  and  all  of  the  moment  derivatives  with  respect 
to  the  same  moment.  We  wish  to  determine  the  stability  characteristics 
of  the  solutions  to  this  set  of  homogeneous,  linear,  constant-coefficient 
differential  equations.  To  do  this  we  assume  solutions  for  each  of 


the  variables  of  the  form  V=V,  (2 
zero  solutions  to  exist  for  the  three  constants 
determinant  of  their  coefficients  must  be  zero. 


etc.  In  order  for  non- 

V  ,  ^  0  ,  the 

Thus, 


(Cxv-vCx^X^ 

X)  (C=5t e  t  0 


=  o 


(193) 
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where  the  lower  subscript  Indicates  the  derivative  with  respect  to  a 
variable  with  the  other  variables  held  constant • 

Expansion  of  this  determinant  gives  a  q,uartlc  eq.uation  in  ^ 

These  four  solutions  for  ^  determine  the  character  of  the  motion 
of  the  aircraft  in  response  to  a  small  disturbance.  We  write  the 
quartic  in  the  form^ 

6  C  0  ^  -V-  E  =  o  (19]+) 

Expansloh  of  the  above  determinant  gives  the  coefficients  of 
this  equation  in  terms  of  the  stability  derivatives.  Eliminating 
the  derivatives  which  are  zero  for 

9-rc  )  —  O  ,  l.e.  for  the  aircraft  Initially  trimmed  level, 
ve  obtain  the  following  expressions  for  the  coefficients: 


A*  Cm©  (195) 

6  =  Cm©  (,C-*Cv 

+  Cm©C-3l^  ~  Cm^  Gx  v  (196) 


C=  Cm©  (Gxv Ciipi  -  -V  Cm©  (G3?.vC;5^ 

^  Qkv  ChiX  (Cz.^Qxq 

-  GxvC^e ^  -  CmcxCxv Ci©  -»■  Cmv  (C^SfCg© 

(197) 
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55“  Ci-iQ  (GscvCz.^-  Cg^Gxc>^)i-  Cm^  Cz.v  Cice 
+  Ch(P<  (^C  gyCxQ  -  Clgy  Czq)  -f  Cmv  (GX(?(C*,q 
~  C^^g(C-XQ  -  Cz^Cx©^ 

E.-  C.3C©  (^C2.vCm(;(~  Cmv/Cx5<^ 


(198) 


(199) 


The  underlined  terms  are  generally  smaller  than  the  others  and  in 
many  cases  can  be  neglected. 

These  expressions  for  the  coefficients  can  be  considerably  simpli¬ 
fied  and  made  more  physically  understandable  by  Introducing  the  con¬ 
cept  of  the  static  and  maneuver  "Msirglns". 

Consider  first  the  airplane  in  trimmed  level  flight.  If  the 
velocity  is  changed  by  a  small  amount,  the  change  in  the  moment  about 
the  Y-axis  with  the  angle  of  attack  free  to  cheuige  but  with  the  re¬ 
striction  that  constant  vertical  force  be  maintained,  i.e.  Cz=  C  ST, 
is  given  by  the  total  derivative, 

cACm  1  _  \  .  ^Cm  \  0(  1  (200) 

dv  'c?.  dv  ^  bo(  dv  IC2L 

Taking  the  derivative  of  the  equation, 

^  (201) 
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we  obtain. 


^Iv 

Solving  this  for  the  total  derivative 


^-o{  I  _,  ^  ^Oit 

dv  Icjt  av 

Thus  we  obtain. 


Tv  Icj,’ 


bc^  Iv 

c^Cm  —  I  _  ^Ch  \ 

dv  bV  diPi 


I, 


or. 


1 

I 

o 

_  -ac,, 

\  dCM  1 

Iv  iCa 

bv 

W  d  C7<  1 

a= 


(202) 


(203) 


(204) 


1  _  ^Ml  ^Cg  fs.H.')  (205) 

IV  ^>v  w  bo^\\/  av  I  I.y  V  7/ 


"5Tiv  ■  av 

This  grouping  of  derivatives  Is  called  the  static  margin,  and 
the  E  coefficient  can  be  expressed  in  terms  of  It  as. 


E=  -  Cxq  (s.m.) 


(206) 


Now  consider  the  aircraft  in  a  steady  pull-up  (See  Fig.  10 ). 

The  change  In  the  moment  about  the  Y-axis  due  to  the  rate  of  rotation 


of  the  aircraft  Is  given  by  the  total  derivative. 

Cm  I  _  b  Om  I  ^Om  I  1 
do  ^0  19  Te'iOz: 


(207) 


From  Figure  10  we  see  that  the  Z-force  in  a  steady  pull-up  is, 

Fzo^-  m  CL  ©  -  mn  V  ©  -  O 


(208) 
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and  its  derivative  with  respect  to  0  is, 


=  —  YT>  V 

^>0 

or  in  non-dimensional  form, 

^  -  -  ^mV 

^6  P  5 

This  gives  a  physical  picture  of  the  derivative  and  is  seen  to 
be  the  first  term  of  the  expression  obtained  previously  for  this 
derivative.  The  other  terms  are  due  to  change  in  the  velocity  at  the 
tall  and  rotor  and  are  usual.ly  smaller  than  the  term  obtained  above. 
For  -  c’st  .  we  obtain, 

do  t  do<=  o 

from  which  we  can  obtain  the  total  derivative, 


(209) 


(210) 


alic/  - 


Thus, 


^0* 


'hc^  Iv 

c^Ch  -  I 

dTS"  'bb  boK  'V 

Multiplying  by  1  /  we  obtain, 

<^Ph  _  'dCiM  _  bQjt  bCH\  = 

bo<  Iv  ^  b(7S  w  ^  Iv 


(M.  M.) 


This  group  of  terms  is  usually  called  the  maneuver  margin  and  appears 
in  the  0  coefficient  and  in  the  ^  coefficient. 


6i 


(211) 


(2J2) 


(213) 


(21k) 


If  ve  make  these  substitutions  in  the  expressions  for  the  coef¬ 


ficients  and  omit  the  terms  which  are  small  compared  to  the  others, 
we  obtain  the  following  expressions: 

,A=  ChqCxvC;!,^  (215) 

I^Hg  Cxv  i-  Ch©Cxv  (216) 

C=  C:3r^  Ch.M.)  -V-  CmgC-XvCs^  (21T) 

P-  Cxv  (m.M.)  - 

■^CMvCCx^Q^e  -  (218) 

t  =  —  CxTq  (219) 


Characteristic  Eq.uatlon  in  Hovering 

When  the  aircraft  is  hovering  the  variable  0(  is  no  longer 
defined  since  the  free  stream  velocity  is  zero.  The  equations  cam 
be  considerably  simplified  due  to  the  fact  that  the  coupling  between 
the  vertical  degree  of  freedom  and  the  other  degrees  of  freedom  is 
very  weak,  i.e.  small  forward  velocities  or  angular  rotations  do  not 
cause  any  appreciable  change  in  the  Z-force  and  small  perturbations 
in  the  vertical  velocity  cause  only  negligibly  small  changes  in  the 
X-force  and  the  moment  about  the  Y-axis.  Thus  the  motion  can  be 
approximated  by  the  following  two  equations: 

(Cxv  -v  Cxv  V  -V-  (Cxe  ^  Cx  ©  )^)  ©  =  O 

Cmv  V  Cm  ©}?)©  = 
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Setting  the  determinant  of  the  coefficients  equal  to  zero  gives  a 
cubic  equation  In  ^  vlth  the  following  coefficients: 


B=  Cm  q 

(221) 

0-  CxvCm©  Cxv  ^ 

(222) 

D=  Oxv  C-Mq  —  Ohv  Cxe 

(223) 

-CmvCxq 

(224) 

Calculation  of  roots  to  the  Quartlc  and  Cubic  Equations  for  Various 
Flight  Configurations 

In  order  to  evaluate  the  stability  for  the  various  flight  con¬ 
figurations  It  is  necessary  to  have  trim  conditions  at  each  tilt 
angle  at  which  the  roots  of  the  characteristic  equation  are  desired. 
These  have  been  evaluated  with  the  aid  of  a  digital  computer  for  the 
Vertol  j6  (Ref.  l^f),  and  eure  shown  plotted  as  a  function  of  q(  and  s/ 
in  Figure  11.  The  velocities  on  these  charts  are  for  a  1:5.2  scale 
model.  The  full  scale  velocity  is  thus  2.28  times  the  velocities  on 
this  scale.  Given  these  conditions,  values  for  the  wing  euid  rotor 
derivatives  are  easily  obtained  by  using  the  charts  In  the  appendix 
and  these  derivatives  allow  calculation  of  the  required  force  and 
moment  derivatives.  These  calculations  were  carried  out  for  hovering 
and  at  Intervals  of  10  degrees  down  to  15  degrees  and  also  for  7 
degrees  which  Is  approximately  the  forward  flight  angle  of  attack. 


63 


(in  forward  flight  the  wing  tilt  angle  is  approximately  eq.ual  to  the 
angle  of  attack. )  The  values  obtained  for  the  derivatives  and  for 
the  coefficients  of  the  stability  q.uartic  are  given  in  the  appendix. 
Using  these  coefficients  the  q.uartics  were  solved  and  the  results 
are  shown  in  Figure  12.  Values  of  the  various  parameters  for  the 
Vertol  76  and  the  lift  curve  slope  assumed  are  given  in  Appendix  B. 

The  values  obtained  for  all  of  the  derivatives  and  variables  at  each 
of  the  tilt  angles  are  given  in  Appendix  D. 

It  should  be  borne  in  mind  that  these  calculations  were  intended 
merely  to  provide  a  qualitative  picture  of  how  the  stability  changes 
while  the  aircraft  moves  through  transition  and  consequently  they  can 
not  be  expected  to  give  accurate  values  for  the  period  and  damping  of 
the  actual  airciaft.  Inspection  of  the  various  contributions  to  the 
coefficients  of  the  quartlc  allows  one  to  see  which  axe  the  main  terns 
and  which  can  usually  be  neglected. 

Approximate  Factorization  of  Quartlc 

Often  much  useful  information  can  be  obtained  by  an  approximate 
factorization  of  the  quartlc.  This  gives  an  expression  for  each  of 
the  roots  in  terms  of  the  stability  derivatives.  If  such  eui  approxi¬ 
mation  is  sufficiently  accurate  and  results  in  expressions  which  are 
not  too  complicated,  it  will  enable  one  to  determine  very  easily  the 
effect  of  any  derivative  on  the  quartlc  roots. 

Since  the  roots  which  represent  the  short-period  oscillation  are 
usually  much  larger  than  the  other  two  for  conventional  airplanes,  we 
assume  that  an  approximation  to  the  short  period  oscillation  can  be 
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obtained  by  dropping  the  last  two  terms  of  the  q.uartic.  Doing  this 
and  retaining  only  the  .largest  te.rms  of  the  first  three  coefficients^ 
we  obtain  the  following  quadratic: 


Ch©  >'^+ Cm©  >  +  =o 

To  obtain  an  approximation  to  the  other  two  roots  we  assume  a 
quadratic  of  the  form  a>?-f  b>  +  e  f  and  multiply  it  times  the 
short-period  approximation.  Then^  equating  coefficients  of  like 
powers  of  ^  ,  gives  five  equations  relating  CL  ,  b  >  C  ^  to 

known  coefficients.  Using  the  firsts  fourth  and  fifth  of  these  we 
obtain  the  following  expressions  for  the  coefficients  of  the  quad¬ 
ratic  irtiich  should  approximate  the  remaining  roots  of  the  quartic: 


C 


Xv 


LcMvC»e-CMeC,,]+Cxv 


Ch.m.) 


(226) 

(227) 

(228) 


In  the  short  period  approximation  the  constant  term  is  directly 
proportional  to  (MM)  whereas  in  the  second  quadratic  both  the  constant 
and  "X  terms  are  Inversely  proportlona,!  to  (MM)  .  We  thus  might 
expect  that  the  large.r  (MM),  the  better  the  approximation. 

To  determine  the  accuracy  of  the  approximation,  these  two  quad¬ 
ratics  were  solved  for  a  number  of  tilt  angles  and  compared  with  so¬ 
lutions  of  the  complete  quartic.  This  comparison  is  illustrated  in 
Figure  13-  We  see  that  for  the  lower  tilt  angles  the  approximation 
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is  fairly  good  and  breaks  down  when  the  roots  become  smaller  as  would 
be  expected.  We  may  also  note  that  the  approximation  did  not  predict 
the  Increased  damping  of  the  7  degree  tilt  angle. 
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DISCUSSION 


Variation  of  quartlc  roots  vlth  tilt  angle 

In  order  to  discuss  VTOL  aircraft  stability  it  is  necessary  to 
know  vhat  types  of  response  are  considered  desirable  and  also  since 
the  desirable  type  can  not  always  be  attained^  we  would  like  to  know 
what  Is  acceptable. 

The  stability  criterion  is  most  conveniently  expressed  by  speci¬ 
fying  the  regions  on  the  complex  plane  into  which  the  roots  of  the 
quartic  which  determine  the  nature  of  the  motion  should  fall.  It 
is  known,  of  course,  that  the  roots  should  fall  into  the  left  half 
plane;  however,  as  is  shown  by  the  hovering  helicopter,  this  condi¬ 
tion  is  not  necessary  in  order  that  an  aircraft  be  capable  of  being 
flown.  An  instability  similar  to  that  ^ich  exists  in  most  hovering 
helicopters  would,  of  course,  never  be  tolerated  in  a  conventional 
airplane.  We  thus  see  that  any  criterion  for  desirable  or  acceptable 
stability  will  depend  on  the  type  of  aircraft  which  is  being  con¬ 
sidered.  It  is  also  true,  as  mentioned  before,  that  handling  quali¬ 
ties  depend  on  more  than  the  values  of  the  complex  frequencies.  The 
amount  of  control  available  to  the  pilot  is  also  very  important. 
Several  criteria  have  recently  been  established  for  VTOL  handling 
qualities  by  means  of  simulators  (See  e.g.  Ref.  I5).  Although  such 
results  will  have  to  be  examined  more  closely  by  actual  flight  test, 
they  do  give  some  indication  of  the  type  of  response  which  is  accepta¬ 
ble  to  the  pilot.  When  the  VTOL  has  developed  considerable  velocity 
in  the  conversion,  then  handling  quality  criteria  established  for 
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conventional  sutsonic  aircraft  can  be  used  as  a  basis  for  establishing 
desirable  VTOL  handling  qualities.  Some  of  the  results  of  these  two 
criteria  are  shown  in  Figure  12. 

In  order  to  obtain  some  indication  of  how  the  stability  charac¬ 
teristics  change  as  the  VTOL  aircraft  goes  through  a  constant  altitude 
conversion,  a  sample  calculation  for  a  typical  aircraft  was  carried 
out  using  the  equations  developed  in  Section  I.  The  results  of  these 
calculations  for  the  Vertol  7^  are  shown  in  Figure  12.  Typical  roots 
for  6U1  airplane  and  hovering  helicopter  are  also  shown.  These  re¬ 
sults  serve  only  as  an  indication  of  how  the  roots  vary  with  tilt 
angle  auid  should  not  be  thought  of  as  giving  an  accurate  picture  of 
the  stability  of  the  Vertol  7^  since  the  effects  of  the  fuselage 
were  omitted  from  the  calculations  and  a  large  number  of  simplifying 
assumptions  were  made  in  order  that  approximate  results  might  be  ob¬ 
tained  without  recourse  to  experimental  data. 

In  the  dlscusssion  which  follows  the  effect  of  the  Important 
stability  derivatives  will  be  considered  as  well  as  the  effect  of 
various  physical  parameters  on  these  derivatives.  V/henever  possible, 
root  locus  plots  are  used  to  show  the  effect  of  Increasing  or  de¬ 
creasing  the  value  of  the  Important  derivatives  starting  from  the 
value  >dilch  was  estimated  for  the  sample  aircraft. 

From  Figure  12  we  see  that  the  hovering  roots  are  quite  similar 
to  those  of  the  typical  hovering  helicopter  shown;  i.e.,  there  is  an 
unstable  oscillation  and  a  real  convergence.  The  oscillation  for  the 
hovering  VTOL  doubles  amplitude  in  about  I.5  seconds  while  the  heli¬ 
copter  takes  3*5  seconds  to  double  amplitude.  According  to  the 
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criteria  of  A'Harrah  and  Kwiatkovski  which  is  shown  in  Figure  12,  the 
VTOL  is  unflyable  in  hovering  while  the  helicopter  roots  shown  fall 
into  the  acceptable  region.  It  will  be  noticed  that  the  damping 
ratios  of  the  two  oscillations  are  nearly  the  same.  It  thus  appears 
that  in  this  regime  damping  ratio  is  not  a  significant  measure  of 
the  handling  qualities. 

The  real  convergence  shown  for  the  VTOL  is  quite  rapid,  having 
a  time  to  half  amplitude  of  only  0.5  seconds.  This  is  seen  to  be 
considerably  more  rapid  than  the  helicopter. 

It  is  interesting  to  see  how  these  hovering  roots  may  be  altered 
by  changing  the  various  derivatives  which  appear  in  the  coefficients 
of  the  cubic  eqtiation  which  determines  the  roots.  The  constemt  coef¬ 
ficient  of  the  hovering  cubic  is, 

-C  fij"  ©  ( 229 ) 

Since  Cxg  Is  &-LwEL^S  !Q6^8t'tl'V6 ^  S^Glr't/jLC  !!LS 

by  the  sign  of  the  derivative  Cmaj-  •  following  root  locus 

sketch  for  variation  of  shows  how  the  stability  characteristics 

of  the  aircraft  change  as  this  derivative  is  changed; 
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We  see  from  this  diagram  that  although  the  time  to  reach  half 
an5)litude  of  the  aperiodic  mode  is  decreased  by  increasing  the  static 
stability,  the  time  to  double  amplitude  of  the  oscillation  is  also 
decreased.  Thus,  making  the  aircraft  statically  more  stable  makes 
it  dynamically  more  unstable.  Decreasing  the  static  stability  has 
just  the  opposite  effect.  The  dyiamic  stability  improves  until  it 
reaches  the  real  axis  \diere  the  roots  combine  to  give  tvo  real  con¬ 
vergences,  one  of  which  soon  becomes  unstable.  The  point  at  which 
it  crosses  the  axis  is  i^ere  O  . 

Probably  the  most  important  derivative  in  the  determination  of 
the  stability  characteristics  of  the  hovering  VTOL  is  the  pitch 
damping  Cm  4  .  The  next  sketch  shows  how  the  cubic  roots  are 

effected  by  variation  of  this  derivative: 


We  see  from  this  that  a  definite  Improvement  in  the  hovering 
roots  can  be  obtained  by  increasing  the  pitch  damping.  The  time 
to  half  amplitude  of  the  real  convergence  is  further  decreased 
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while  at  the  same  time  the  damping  characteristics  of  the  oscillation 
are  improved.  This  derivative  in  Ijovering  is  mainly  due  to  the  offset 

4 

hinges  of  the  rotor  and  can  be  Increased  by  increasing  the  eunount  of 
offset . 

The  effect  of  changing  the  moment  of  inertia  about  the  Y-aocls  is 
shown  below: 


Increasing  the  moment  of  inertia  causes  all  of  the  roots  to  move 
toward  the  origin.  By  comparing  this  sketch  with  Figure  12  we  see 
that  Increasing  the  moment  of  inertia  of  the  VTOL  a  sufficient  amount 
makes  them  quite  similar  to  those  of  the  hovering  helicopter  shown  in 
the  drawing. 

We  now  consider  -vdiat  happens  as  the  wing  is  tilted  forward  from 
the  hovering  position.  The  first  fact  which  can  be  noticed  from  an 
examination  of  Figure  12  is  that  the  stability  characteristics  change 
quite  rapidly  as  the  aircraft  begins  to  move  forward.  This  is  due  to 
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the  fact  that  as  soon  as  the  aircraft  has  some  forward  velocity  the 
tail  has  an  important  effect  on  the  derivatives  whereas  in  hovering 
it  had  practically  no  effect  at  all. 

With  ten  degrees  of  tilt  angle  the  short  period  oscillation 
changes  period  only  slightly,  but  the  damping  characteristics  are 
Improved  considerably.  As  the  wing  is  tilted  down  another  10  degrees 
the  real  part  of  the  root  changes  from  positive  to  negative  while 
the  period  decreases  to  about  four  seconds. 

It  is  worth  mentioning  here  that  there  is  no  definite  reason 
why  the  loci  sketched  through  the  points  in  Figure  12  should  be 
drawn  as  they  are.  Thus,  it  is  possible  that  the  mode  which  becomes 
similar  to  the  airplane  short  period  could  be  associated  with  the 
mode  which  is  represented  by  the  two  real  convergences  near  hovering. 
Thus,  the  locus  through  the  roots  could  appear  as  in  the  following 
sketch : 
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The  dotted  lines  show  where  the  locus  was  drawn  in  Figure  12.  Which 
of  these  lines  is  correct  could  be  determined  by  computing  the  roots 
for  several  more  values  of  the  tilt  angle  between  65  and  75  degrees. 
As,  however,  this  has  not  yet  been  done,  we  will  assume  in  the  follow¬ 
ing  discussion  that  the  locus  is  as  shown  in  Figure  12. 

When  the  aircraft  attains  some  forward  velocity  the  Z-force  e- 
quation  is  no  longer  uncoupled  from  the  X-force  and  pitching  moment 
equations  and  there  are  then  fotir  roots  rather  than  three.  The 
fourth  root  was,  of  course,  present  in  hovering  but  it  was  not  neces¬ 
sary  to  compute  it  since  it  is  uncoupled  from  the  others  and  is  a 
very  large  negative  root  and  thus  has  little  effect  on  the  motion  of 
the  aircraft. 

As  the  aircraft  gains  some  velocity  this  root  moves  in  along  the 
real  6txls  and,  as  shown  by  Figure  12,  combines  with  the  other  real 
convergent  root  to  form  a  long  period  damped  oscillation  which  be¬ 
comes  more  lightly  damped  as  the  wing  tilts  further  forward.  Between 
70  and  65  degrees  this  oscillation  again  becomes  two  real  roots,  one 
of  which  becomes  unstable.  We  thus  see  that  a  rather  radical  change 
tedtes  place  in  both  of  the  characteristic  modes  at  a  tilt  angle  of 
about  65  degrees.  This  is  due  to  the  fact  that  the  wing  stalls  at 
this  point  in  the  conversion. 

At  this  point  it  is  useful  to  see  what  happened  to  the  equation 
to  cause  the  stability  characteristics  to  deteriorate  so  badly.  To  do 
this  we  consider  the  roots  for  a  tilt  angle  of  65  degrees  given  in 
Figure  12  and  use  the  root  locus  plot  to  determine  what  changes  in  the 
static  and  maneuver  margin  would  be  necessary  to  obtain  satisfactory 
or  at  least  improved  stability  characteristics  at  this  tilt  angle. 
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For  this  tilt  angle  we  have  a  dynamically  stable  oscillation 
and  two  real  roots,  one  of  which  decreases  to  half  amplitude  in  2.8 
seconds  while  the  other  doubles  amplitude  in  2.3  seconds.  In  order 
to  obtain  satisfactory  flying  qualities,  the  divergent  root  should 
be  removed  or  at  least  moved  much  closer  to  the  origin,  i.e.  the  time 
required  to  double  amplitude  should  be  increased.  As  the  damping  is 
only  marginal,  it  should  be  increased. 

In  Figure  ih  we  see  the  effect  of  changing  the  static  margin. 

The  gains  shown  on  this  locus  are  not  for  changes  in  the  static 
margin  but  for  -  times  the  static  margin;  where  A  is 

the  coefficient  of  the  term  in  the  quartic  equation  in  ^  . 

If  (SM)  is  decreased  to  zero  the  two  real  roots  come  together.  For 
further  decrease  in  (SM)  the  roots  leave  the  real  axis  to  form  a 
sHghtty- unstable  oscillation  with  a  time  to  double  amplitude  of  about 
17  seconds.  For  this  amount  of  change  in  (SM)  we  see  that  the  short 
period  oscillation  is  hardly  effected  at  all.  Further  decrease  in 
the  static  stability  causes  the  period  of  both  oscillations  to  de¬ 
crease.  One  of  the  oscillations  becomes  more  unstable  while  the 
other  becomes  more  stable.  We  thus  see  that  making  the  aircraft 
statically  stable  improves  its  characteristics  but  that  too  much 
static  stability  is  undesirable  since  it  results  in  a  very  short- 
period  divergent  oscillation.  We  may  also  note  that  no  amount  of 
change  in  the  static  stability  alone  will  make  the  stability  charac¬ 
teristics  completely  satisfactory,  i.e.  make  both  modes  both  stati¬ 
cally  and  dynamically  stable  with  a  desirable  damping  ratio. 
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In  Figure  I5  is  shovn  the  result  of  changing  (MM).  As  (MM)  is 
increased  the  real  divergence  moves  closer  to  the  origin;  however^ 
as  may  he  seen  from  the  values  of  the  gain  shown,  large  gain  changes 
result  in  only  a  small  movement  while  the  period  of ' the  short -period 
mode  decreases  rapidly  with  Increase  in  the  gain.  Since  the  real 
part  of  the  oscillatory  root  hardly  changes,  the  damping  ratio  is  de¬ 
creased.  Decreasing  (MM),  while  it  improves  the  damping  ratio  of  the 
short -period  mode,  decreases  the  time  to  double  amplitude  of  the  real 
divergence  and  thus  is  not  a  desirable  change.  We  thus  see  that 
little  if  any  Improvement  can  be  obtained  by  changing  (MM). 

Figure  16  shows  the  effect  of  changing  ,  which  appears 

in  both  (SM)  and  (MM).  Decreasing  Ohlo<  (increasing  1  OmcK  I  ) 
decreases  the  period  of  the  short-period  mode  and  has  little  effect 
on  the  real  roots .  We  see  however  that  Increasing  Ofn(A  iJuproves 
the  stability  characteristics  considerably.  If  it  is  increased 
enough  to  make  the  aircraft  statically  stable,  then  there  will  be  a 
long-period  oscillation  with  a  time  to  double  amplitude  of  about  10 
seconds.  This  amount  of  change  in  Cl(m(^  also  improves  the  short 
period  oscillation  since  the  damping  ratio  is  increased  considerably. 
According  to  this  root  locus,  if  the  value  of  were  increased 

from  -.439  to  about  -.14  the  stability  characteristics  would  be  con¬ 
siderably  Improved.  It  must  be  realized,  however,  that  changing 

in  this  manner  Improves  the  stability  only  because  the  wing 
is  stalled  at  this  tilt  angle.  The  derivative  appears  in  the 

static,  margin  multiplied  by  Qzv/  which  is  almost  always  negative 
when  the  wing  is  not  stalled.  When  the  wing  is  stalled  it  can  be 
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positive  as  it  is  in  this  case.  This  makes  the  effect  on  the  (SM) 
of  changing  Cnr>(^  just  the  opposite  of  what  it  would  he  for  an 
unstable  wing.  For  all  considerations  such  as  these  it  should  he 
remembered  that  the  conditions  at  which  the  wing  will  stall  can  not 
he  predicted  accurately  in  the  power-on  case.  The  wing  was  assumed 
to  obey  the  power  off  lift  curve  shown  in  Figure  39  a  number  of 
effects  are  present  which  may  cause^the  stall  to  occur  at  a  different 
effective  angle  of  attack. 

As  the  wing  is  tilted  further  forward  from  65  degrees^  the 
effective  angle  of  attack  becomes  still  greater  and  the  wing  remains 
staJ-led.  Figure  I8  shows  how  the  effective  angle  of  attack  varies 
through  the  conversion.  This  chart  may  be  used  to  quickly  find  the 
effective  angle  of  attack  at  any  tilt  angle  if  the  forward  velocity 
and  the  thrust  coefficient  are  known.  As  may  be  seen  by  this  chart 
the  angle  of  attack  increases  with  increase  in  the  velocity  and  de¬ 
creases  with  Increase  in  the  thrust  coefficient.  As  a  consequence 
of  this  fact,  the  effective  angles  of  attack  will  be  lower  for  an 
accelerating  aircraft  and  higher  for  a  decelerating  aircraft.  The 
sample  conversion  shown  on  this  chart  is  for  equilibrium  trim  con¬ 
ditions,  i.e.  the  trim  conditions  were  computed  at  each  tilt  angle 
as  if  the  airplane  were  in  equllibrltim  flight  at  this  tilt  angle. 

In  Figure  I9  is  shown  a  plot  of  effective  angle  of  attack  versus 
tilt  euagle  for  different  amounts  of  acceleration  and  deceleration. 

We  see  from  this  that  even  small  amounts  of  acceleration  and  deceler¬ 
ation  can  change  the  angle  of  attack  considerably.  Also  note  that 
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the  tilt  angle  at  which  the  maximum  angle  of  attack  is  attained  in¬ 
creases  for  a  decelerating  airciaft  and  decreases  for  an  accelerating 
aircraft.  We  also  see  that  for  this  particular  aircraft^  an  acceler¬ 
ation  of  about  0.2  g  is  required  in  order  for  the  wing  to  he  unstalled 
over  the  entire  conversion. 

The  lift  coefficients  obtained  using  angles  of  attack  obtained 
in  this  way  and  the  power  off  lift  curve  slope  should  not  be  expected 
to  be  too  accurate  due  to  the  assumptions  which  have  been  made.  The 
following  three  effects  have  been  found  experimentally  to  be  present 
in  the  slipstream  (Ref.  4)  and  may  cause  the  lift  developed  to  differ 
somewhat  from  that  predicted  on  the  basis  of  simplifying  assumptions: 

1 .  Slipstream  rotation . 

2.  Delayed  stall  due  to  an  effective  boundary  layer  control. 

3.  Variation  in  induced  velocity  across  disc, 

Brenkmann  concludes  that  for  prediction  of  the  average  lift  on  the 
wing  the  delayed  stall  effect  is  the  most  importemt  and  that  potential 
flow  solutions  may  give  satisfactory  results  if  this  additional  lift 
is  taken  into  account.  He  also  shows  that  for  the  wing -propeller 
combination  which  he  tested,  slender -body  theory  (See  Ref.  5)  predicts 
the  lift  very  well  if  the  contribution  due  to  the  de.layed  stall  is 
added  to  that  given  by  the  theory. 

As  the  tilt  angle  decreases  further  we  see  that  the  period  of 
the  short -period  mode  decreases  with  essentially  no  change  in  the 
time  of  half  eunplitude.  This  means,  however,  that  the  damping  ratio 
of  the  mode  is  decreased.  The  period  of  this  mode  thus  changes  from 
4.2  seconds  at  a  tilt  angle  of  65  degrees  to  2,2  seconds  at  a  tilt 


77 


angle  of  15  degrees.  According  to  the  boundaries  established  by 
means  of  an  F-9^,  (Ref.  l6)  the  damping  of  this  mode  is  still  not 
sufficient,  although  it  is  not  too  far  from  the  acceptable  boundary. 

With  further  decrease  in  the  tilt  angle  the  other  mode  comes 
back  to  the  real  axis,  giving  for  a  tilt  angle  of  65  degrees  a  pure 
divergence.  With  further  decrease  the  time  to  double  amplitude  of 
the  divergent  mode  at  first  becomes  shorter  and  then  longer  again. 

We  see  that  these  roots  remain  on  the  real  axis  essentially  until 
the  angle  of  attack  becomes  small  enough  for  the  ving  to  become  un¬ 
stalled. 

We  now  consider  how  the  roots  in  this  region  will  be  effected 
by  variation  in  (SM)  and  (MM)  as  we  did  before  for  a  tilt  angle  of 
65  degrees.  In  this  case  as  a  representative  angle  we  will  use  a 
tilt  angle  of  35  degrees. 

In  Figure  17  we  see  the  results  of  changing  (SM)  starting  from 
a  tilt  angle  of  35  degrees.  As  it  is  decreased  the  two  real  roots 
move  together,  the  divergent  root  meeting  the  imaginary  axis  when 
(SM)  becomes  zero.  As  (SM)  is  further  decreased  the  two  roots  leave 
the  real  axis  to  form  a  stable  oscillation.  We  note  that  although 
the  roots  to  begin  with  are  quite  similar  to  those  in  the  65  degree 
tilt  angle  case,  the  oscillation  here  is  convergent  whereas  in  the 
other  case  it  was  slightly  unstable.  For  gain  sufficient  to  keep 
this  mode  from  crossing  back  over  the  Imaginary  axis  the  short-period 
mode  has  moved  only  a  small  distance  and  in  the  direction  of  increased 
damping.  We  then  see  that  for  this  tilt  angle,  giving  the  aircraft  an 
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adequate  amount  of  static  stability  makes  the  period  and  damping  of 
both  modes  acceptable. 

Changing  (MM)  for  this  tilt  angle  is  almost  the  same  as  for  the 
65  degree  tilt  angle  and  is  thus  not  shown. 

o 

In  Figure  20  we  see  the  effect  of  changing  ,  the  deriva- 

tiye  of  the  pitching  moment  with  respect  to  angle  of  attack  for  constant 
velocity.  We  see  that  decreasing  Ohct^  almost  no  effect  on  the 

two  real  roots  since  the  zeros  for  this  root  locus  are  located  almost 
bn  these  poles.  The  period  of  the  short-period  mode  is  decreased  with 
etlmost  no  change  in  time  to  half  amplitude  for  a  decrease  in  Cndi  » 
Increasing  Oh^A  some'vrtiat  more  interesting.  The  real  divergent 
root  moves  toward  the  Imaginary  axis  and  the  real  convergent  root  de¬ 
creases  its  time  to  half  amplitude.  We  see  however  by  comparing  the 
gains  on  the  two  loci  that  a  Ch  d  increase  sufficient  to  make  the 
divergent  mode  stable  causes  a  very  large  movement  of  the  short-period 
mode.  In  fact  at  about  the  same  gain  where  the  one  mode  makes  the  air¬ 
craft  statically  stable,  the  other  mode  becomes  dynamically  unstable. 
Thus,  changing  CTj-y  cA  improve  the  dynamic  characteristics  for  this 
tilt  angle  somewhat  although  it  will  not  make  them  satisfactory. 

Finally,  going  back  to  Figure  12,  we  see  what  happens  to  the  two 
modes  when  the  tilt  angle  becomes  small  enough  to  unstall  the  wing. 

In  this  case  the  tilt  angle  is  approximately  equal  to  the  effective 
angle  of  attack.  The  angle  of  attack  in  normal  forward  flight  is  about 
7  degrees  and  this  tilt  angle  is  shown  in  the  diagram.  The  damping  of 
the  short  period  mode  decreases  conslderab.ly  without  much  change  in 
the  period  and  the  two  real  roots  form  a  long  period  oscillation 


79 


similar  to  the  usual  phugold.  The  aircraft  in  forward  flight  then  has 
acceptable,  although  not  ideal,  stability  characteristics. 

Effect  of  Riysical  Rsuumeters 

In  the  last  section  we  saw  how  the  q.uartic  roots  varied  as  the 
tilt  angle  was  changed  for  a  typical  tilt -wing  VTOL  aircraft.  It  was 
also  shown  by  means  of  root  locus  plots  how  the  stability  character¬ 
istics  at  several  tilt  angles  would  be  changed  if  the  values  of  the 
two  groups  of  terms,  (SM)  and  (MM),  were  changed.  We  would  now  like 
to  see  how  the  various  derivatives  >diich  make  up  these  two  groups  de¬ 
pend  on  the  aircraft  physical  parameters. 

In  our  original  perturbatlcn  equations,  there  were  l8  deriva¬ 
tives.  Six  of  these  were  found  to  be  either  zero  or  negligibly  small. 
Of  the  remaining  derivatives,  five, 

,  depend  essentially  on  only  the  mass,  moment  of  inertia,  and 
forward  velocity.  Two  of  these  derivatives,  C^Q  >  ,  also 

have  small  contributions  from  the  tail  and  other  parts  of  the  aircraft, 
but  these  are  small  compared  to  the  inertia  contributions.  Thus  we 
will  be  mainly  Interested  in  the  derivatives  of  the  forces  and  moment 
with  respect  to  eoigle  of  attack  and  free  stream  velocity  and  the  pitch 
damping  derivative  0  q  .  All  of  these  seven  derivatives  depend 
rather  strongly  on  the  easily  alterable  physical  parameters  of  the  air¬ 
craft  and  are  thus  of  greatest  Interest. 

Consider  first  the  pitch  damping.  This  derivative  is  mainly  due 
to  the  tall  although,  as  will  be  seen,  the  rotor  can  also  be  quite 
•  Important.  The  fuselage  also  contributes  some  pitch  damping  but  this 
can  not  be  expressed  in  analytical  form  for  fuselages  of  unusual  shape 
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and  must  therefore  he  determined  experimentally  for  a  particu.lar  con¬ 
figuration.  Also,  the  pitch  damping  of  the  fuselage  is  usually  quite 
small  compared  to  that  from  the  other  parts  of  the  aircraft. 

The  contributions  to  the  pitch  damping  as  a  function  of  the  rotor 
angle  of  attack  are  shown  in  Figure  21.  We  see  that  the  largest  con¬ 
tribution  here  is  that  of  the  horizontal  stabilizer.  This  contri¬ 
bution  is  simply  a  result  of  the  fact  that  when  the  aircraft  is  pitch- 
Ing  at  a  rate  ©  ,  the  angle  of  attack  ef  the  tail  is  increased  by  an 

ie 

amount  ,  producing  a  moment  which  opposes  the  pitch  rate.  This 

moment  is  thus  directly  proportional  to  the  distance  from  the  center 
of  gravity  to  the  tall  and  to  the  velocity  of  the  fluid  at  the  tall. 
The  contribution  of  the  offset  hinges  is  seen  to  be  independent  of 
the  angle  of  attack  and  thus  of  the  free  stream  velocity.  It  is  for 
this  reason  that  this  contribution  is  particularly  important.  Al¬ 
though  it  is  small  compared  to  the  tail  contribution  in  forward  flight, 
it  is  almost  the  only  source  of  damping  in  hovering.  This  contri¬ 
bution  is  directly  proportionai  to  the  amount  of  hinge  offset.  We 
can  see  how  this  moment  arises  in  the  following  way: 

Consider  the  side  view  of  a  two  bladed  rotor  with  a  mean  coning 
angle  |So  «uid  a  flapping  angle  0-i  (See  Fig.  7), 
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The  moments  about  points  A  &  B  due  to  the  lift  on  the  rotor  and  the 
centrifugal  force  must  balance  since  the  blades  are  hinged  at  these 


points .  Thus , 

%  ^  z  ~~  X  CLP  -  O 


(230) 


and 


Fz  =  CF  "tojn  P  =  CF  p  (231) 

since  p  Is  small.  The  force  F ^  produces  a  moment  about  point 

0  due  to  the  offset  e.  The  moment  due  to  the  two  blades  Is  thus^ 

K\=  -CFe[^Cp>o-G^0-(Bo-v-o.,)1[-=  aCF  <2  (232) 

Thus,  a  moment  Is  produced  which  Is  proportional  to  C3L»j  .  When  the 
aircraft  Is  rotated  at  a  rate  ©  about  Its  center  of  gravity  the 
rotor  can  not  retain  Its  former  position  since  It  Is  hinged  and  can 
only  be  moved  by  aerodynamic  forces.  The  rotor  thus  lags  behind 
producing  em  effective  change  In  O-j  of  magnitude. 


bOL, 


U  ) 


(233) 


We  thus  have, 

^  S>.CF  C-'-Sfe)  (234) 

o©  0^ 

In  the  actual  case  the  rotors  are  moving  and  these  forces  must  be 
averaged  over  the  rotation  angle  'F  was  done  previously.  The 
simple  considerations  above  show,  however,  the  source  of  the  velocity 


independent  moment. 
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The  contribution  of  the  tail  fan  although  non-flapping  l,s  also 
independent  of  the  velocity  and  is  very  small  due  to  the  smallness 
of  this  rotor. 

The  contribution  of  the  rotor  forces  are  particularly  inter¬ 
esting  due  to  the  fact  that  above  a  certain  velocity  the  H-force 
actually  contributes  a  "negative  damping"  force.  This  phenomenon 
can  be  explained  in  the  following  way:  When  the  aircraft  is  pitching 
at  a  given  rate,  the  rotor  lags  behind  its  position  for  the  non¬ 
pitching  rotor  as  mentioned  previously.  Due  to  this  lage  the  blades 
on  one  side  of  the  rotor  are  flapping  up  and  on  the  other  side  are 
flapping  down.  The  velocity  components  due  to  this  flapping  causes 
the  angle  of  attack  of  the  b.lades  to  differ  on  the  two  sides  and 
hence  to  change  the  magnitude  and  direction  of  the  local  lift  vector. 
When  the  Inflow  is  small,  the  change  in  direction  is  greater  than  the 
magnitude  change  and  results  in  an  H-force  change  which  opposes  the 
pitch  rate.  However,  when  the  inflow  becomes  sufficiently  large, 
the  magnitude  change,  which  is  in  the  opposite  direction,  is  the 
greatest  and  the  H-force  change  tends  to  aid  the  pitch  rate.  The 
pitch  damping  due  to  thrust  change  is  then  positive.  Tliis  effect  is 
illustrated  in  Figure  22. 

We  see  that  the  pitch  damping  can  be  increased  by  increasing  the 
tail  size  and  the  amount  of  flapping  hinge  offset,  and  that  there 
exists  a  destabilizing  contribution  to  the  pitch  danQ)ing  which  in¬ 
creases  with  the  inflow  ratio. 

Consider  next  the  moment  resulting  from  an  angle  of  attack  change 
at  constant  velocity.  The  important  contributions  to  this  derivative 
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are  shown  In  Figure  23.  We  see  that  the  largest  contribution  here 
is  definitely  that  due  to  the  change  in  angle  of  attack  of  the  tail 
and  it  depends  on  the  distance  between  the  aircraft  center  of  gravity 
and  the  tail  and  on  the  fluid  velocity  at  the  tail.  All  of  the  other 
contributions  are  smll  but  desta^bilizing  except  that  of  the  wing 
Z-force.  This,  however,  is  due  to  the  fact  that  the  wing  is  stalled. 

If  the  wing  is  not  stalled  throughout  the  transition,  all  contri¬ 
butions  to  this  derivative  except  the  tail  are  always  destabilizing. 
When  the  wing  is  not  stalled  there  is  a  center  of  gravity  position 
at  which  the  wing  and  rotor  contribution  cancel  that  of  the  tail. 

The  center  of  gravity  range  is  thus  decreased  due  to  the  destabili¬ 
zing  contributions  of  the  rotor. 

We  must  remember,  however,  that  positive  does  not  neces¬ 

sarily  mean  static  instability  as  is  usually  the  case  with  conventional 
airplanes.  For  the  VTOL  the  static  stability  is  determined  by  the 
sign  of, 

(SM)=  Cx«C  M\/  ”  Czv/  (235) 

and  the  first  term  may  be  of  equal  or  greater  importance  than  the 
second.  The  relative  magnitude  of  these  two  terms  as  a  function  of 
rotor  auigle  of  attack  is  shown  in  Figtire  2k.  From  this  plot  we  see 
that  near  hovering  both  terms  are  stabilizing.  At  about  20  degrees 
the  wing  stalls  and  both  terms  change  sign  and  become  destabili¬ 
zing.  The  C^y/  term  is  then  very  small  up  to  about  80  degrees 
while  the  Cm  term  is  destabilizing  in  this  region.  At  about 
83  degrees  both  terms  again  become  stabilizing.  Thus,  we  see  that 
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the  derivative  ''-MV  is  quite  important  near  hovering  and  in 
forward  flight. 

The  importsuit  contributions  to  CM^/  are  shown  in  Figure  25. 

The  major  contribution  in  hovering  is  the  wing  force  in  the  X-direction. 
This  derivative  is  mainly  due  to  the  change  in  the  lift  of  the  wing 
in  the  slipstream  since  the  wing  is  in  a  nearly  vertical  position. 

It  owes  its  magnitude  to  the  fact  that  there  is  considerable  verti¬ 
cal  separation  between  the  wing  center  of  pressure  and  the  center  of 
gravity  of  the  aircraft,  which  is  not  present  in  a  conventional  air¬ 
plane.  When  the  single  of  attack  is  high  enough,  however,  the  wing 
stalls  and  this  contribution  becomes  small  and  of  opposite  sign.  The 
two  rotoi'  terms  we  see  depend  on  the  relative  position  of  the  rotor 
and  the  aircraft  center  of  gravity.  The  overall  effect  of  the  rotor 
is  to  give  a  small  positive  contribution  to  Om\/  which  increases 
in  magnitude  as  the  wing  is  tilted  down. 

Positive  Omv/  usually  Increases  the  static  stability  since 
is  always  negative  for  a  non-stalled  wing  and  (SM)  is  nega¬ 
tive  for  static  stability.  The  following  shows  how  the  signs  of  the 
derivatives  effect  (SM), 


(5^)  ~  Cmv 


c 


M  o( 


(  -  )  = 

Usual  sign 
for  unstalled, 
wing. 


May  change  sign  for  stalled 
wing. 


(  -  ) 

Always  negative 
if  tail  is  suf¬ 
ficiently  large. 
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We  see  from  this  that  if  wing  stall  causes  to  change 

sign,  this  wl.ll  give  a  destabilizing  contribution.  If  both 
and  Omv/  change  the  contribution  will  still  be  stabilizing  while 
if  either  changes  sign  while  the  other  does  not,  the  contribution 
will  be  destabilizing. 

In  Figure  26  we  see  the  contributions  to  ,  the  change 

in  force  in  the  Z-direction  with  angle  of  attack  at  constant  velocity. 
In  a  conventional  aircraft,  if  the  tail  lift  is  neglected,  this  deriva¬ 
tive  is  just  the  slope  of  the  wing  lift  curve.  Here  we  see  that  while 
the  effect  of  the  wing  is  still  Important,  the  rotor  also  gives  a 
large  contribution.  Although  the  stalled  wing  gives  a  large  positive 
contribution,  the  negative  contributions  of  the  rotor  and  the  tail 
are  almost  stifficlent  to  keep  this  derivative  negative  throughout 
conversion.  The  tall  contribution  is  simply  due  to  the  change  in 
tall  lift  wl.th  angle  of  attack  and  hence  is  proportional  to  the  dy¬ 
namic  pressure  at  the  tail.  The  triangular  shape  of  the  rotor  contri¬ 
bution  can  be  explained  as  follows:  If  the  angle  of  attack  is 
slightly  changed,  this  effects  the  Inflow  to  the  rotor,  ^  ,  £ind 

hence  the  thrust.  Thus,  near  hovering  there  is  little  contribution 
to  since  the  rotor  is  nearly  parallel  to  the  flow  and  the 

inflow  is  hardly  effected  by  a  small  angle  of  attank  change.  In  for¬ 
ward  flight,  on  the  other  hand,  the  thrust  vector  is  almost  parallel 
to  the  flow  and  the  thrust  force  is  in  the  X-direction  and  does  not 
effect  .  Since  both  these  effects  vary  approximately  linearly 

with  tilt  angle,  the  rotor  will  have  a  maximum  Z-force  change  with 
angle  of  attack  at  a  tilt  angle  of  about  45  degrees. 


86 


the 


The  other  derivative  which  appears  in(SM)is  > 

change  in  Z-force  with  velocity  at  constant  angle  of  attack.  This 
derivative  has  a  contribution  from  the  rotor  and  from,  the  wing. 

These  are  shown  in  Figure  27.  As  mentioned  before,  this  derivative 
is  usustUy  negative  for  an  unstalled  wing.  However,  we  see  that  the 
rotor  contribution  is  positive  euid  fairly  large  and  hence  could  be 
large  enough  to  make  a  v/  positive  even  for  an  unstall.ed  wing. 

The  .rotor  tem  arises  from  the  fact  that  if  the  velocity  in¬ 
creases  at  constant  angle  of  attack  the  inflow  to  the  rotor  is  in¬ 
creased,  decreasing  the  angle  of  attack  of  the  rotor  b.lades  and  de¬ 
creasing  the  thrust.  The  thanist  decrease  becomes  greater  as  the  wing 
is  tilted  down  and  the  freestream  velocity  has  a  greater  component 
perpendicular  to  the  rotor  disc.  Near  the  forward  fll^t  position 
the  thrust  is  mostly  in  the  X-direction  and  the  contribution  to  CWv/ 
is  decreased. 

The  wing  contribution  here  depends  on  wing  angle  of  attack  and 
becomes  positive  or  destabilizing  when  the  wi.ng  stalls. 

Consider  now  how  the  derivatives  we  have  discussed  effect  the 
maneuver  margin.  The  signs  of  the  various  contributions  are  shown 
by  the .following; 


(  +  ) 

Usual  sign 
for  unstalled 
wing. 


(  -  ) 

May  change  Always  same  sign, 
sign  for 
stalled  wing. 


87 


We  see  from  this  that  the  wing  stall  can  make  the  first  term  decrease 
the  value  of  (MM) .  An  increase  in  the  magnitude  of  all  the  other 
terms  tends  to  increase  (MM).  The  importance  of  (MM)  cam  be  seen 
from  the  short  period  approximation  given  earlier. 


(236) 


where  (MM)  in  this  expression  is  not  made  non-dimensional.  From 
this  we  see  that  Increasing  (MM),  decreases  the  period  of  the  short 
period  mode  am  was  also  shown  by  the  root  locus  plots  in  the  previous 
section. 

We  will  consider  now  the  two  X-force  derivatives.  These  deriva¬ 
tives  do  not  appear  in  either  (SM)  or  (MM)  and  do  not  effect  the 
stability  characteristics  as  much  as  those  derivatives  already  dis¬ 
cussed. 

The  change  in  the  force  in  the  X-directlon  with  an  increase  in 
angle  of  attack  is  normally  negative,  i.e.  if  the  angle  of  attack  is 
increased,  a  force  is  produced  which  tends  to  oppose  the  motion  of  the 
aircraft.  The  contributions  to  this  derivative  are  shown  in  Figure  28. 
The  tall  contribution,  as  wotild  be  expected,  is  small,  negative  and 
increases  in  proportion  to  the  velocity.  The  rotor  contribution  is 
positive  and  is  given  approximately  by 

Thus,  an  increase  in  angle  of  attack  causes  an  increased  inflow  which 
increases  the  thrust,  producing  a  positive  force  in  the  X-dlrection. 

The  effect  of  the  thrust  Increase  becomes  greater  as  the  wing  is 
tilted  forward.  The  various  terms  in  the  wing  contribution  are  shown 
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in  the  upper  part  of  this  Figure.  We  see  that  the  biggest  term 

is  that  involving  a  change  in  the  angle  of  attack,  i.e.  (  \ 

rather  than  a  change  in  the  magnitude  of  the  resultant  velocity. 

The  fact  that  the  wing  is  stalled  decreases  the  value  of  this  contri¬ 
bution  but  the  term  involving  the  lift  coefficient  is  larger  at  small 

tilt  angles,  and  the  stalled  wing  has  only  a  small  effect  on  the  con- 

tributioi;.  The  sum  of  the  various  contributions  is  negative  as  in  a 
normal  aircraft  and  its  magnitude  Increases  as  the  tilt  angle  decreases. 

The  X-force  change  with  velocity  is  also  normally  negative  since 
a  velocity  increase  will  usually  produce  a  force  tending  to  oppose 
the  motion  of  the  aircraft.  The  rotor  contribution  to  the  derivative 

is  always  negative  and  approximately  proportional  to  . 

O/Mx 

is  always  negative,  i.e.  an  increase  in  the  freestream  ve¬ 
locity  produces  a  decrease  in  the  inflow  to  the  rotor.  This  causes 
a  decrease  in  the  thrust  and  hence  a  force  in  the  negative  X-directlon. 
The  effect  of  this  thrust  decrease  increases  as  the  tilt  angle  is  de¬ 
creased.  This  is  shown  in  Figure  29. 

The  various  terms  in  the  wing  contribution  are  also  shown  in  this 
Figure.  Again  the  main  effect  is  the  change  in  angle  of  attack  rather 
than  the  change  in  the  magnitude  of  the  resultant  velocity.  The  sharp 
drop  at  about  70  degrees  is  due  to  wing  stall.  The  sum  of  the  wing 
and  rotor  contributions  has  a  dip  at  about  70  degrees  since  the  nega¬ 
tive  rotor  term  is  still  small  at  this  tilt  angle. 
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coNcnjsioiK 


1.  Values  for  the  deriratiyes  of  the  resultant  velocity  and 
effective  angle  of  attack  with  respect  to  forward  velocity  and  air¬ 
craft  angle  of  attack  can  he  obtained  from  general  charts  without 
extensive  calculations . 

2.  Typical  tilt  wing  VTOL  has  a  short  period  mode  whose  period 
becomes  smaller  as  the  tilt  angle  is  decreased.  The  time  to  half 
amplitude  of  this  mode  is  approximately  equal  to  a  constant  times 
the  pitching  moment  of  inertia  divided  by  the  pitch  damping.  The 
other  mode  is  a  long  period  oscillation  which  becomes  a  pure  con¬ 
vergence  and  a  pure  divergence  when  the  wing  stalls. 

3.  For  unaccelerated  level  flight  the  wing  stalls  over  part 

of  the  transition.  The  amount  of  stall  is  decreased  by  acceleration 
or  climb  and  increases  by  deceleration  or  descent. 

4.  The  aircraft  very  quickly  becomes  statically  unstable  when 
the  wing  stalls. 


RBCOIttaaDATIOMS  FOR  FURTHER  STUDY 

It.  order  to  gain  a  better  xmderstanding  of  VTOL  dynamics,  it  is 
necessary  that  more  experlaeatal  work  be  done.  Detailed  analytical 
analysis  are  not  particularly  useful  if  it  cannot  be  detennlned 
whether  or  not  they  accurately  represent  the  physical  situation.  Two 
types  of  experimental  work  are  needed.  First,  it  would  be  desirable 
to  have  accurate  measurements  of  the  forces  on  the  wing  due  to  the 
slipstream  in  order  to  determine  idiat  parameters  of  a  rotor  wing 


combination  effect  the  lift  oiwracterlstics .  Also,  d;ynamic  stability 
information  is  needed  to  detemlne  vhether  oi  not  present  theoretical 
treatments  give  an  accurate  description  of  VTOL  dynamic  characteristlcc . 

A  theoretical  treatnaEtt  of  the  dynamic*  of  a  complete  tremsition 
vith  time-varying  stability  derivatives  voxild  also  be  valuable  since 
this  type  of  analysis  nd^t  differ  considerably  from  the  steady  state 
analysis  If  the  acceleration  is  aafficiently  high. 
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Figure  28 
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APPENDICES 


A.  Inflow  Ratio  and  Induced  Velocity  Derivatives 

In  this  section  we  will  give  expressions  for  the  derivatives  of 
y\  ,  the  inflow  ratio,  and  KT  ,  the  induced  velocity,  with  respect 
■to  CK  ,  fJi'^  )  and  .  These  expressions  are  also  plotted  in 

the  form  of  charts  in  Figures  30  through  37-  The  H-force  coefficient 
derivatives  which  were  derived  earlier  are  plotted  in  Figure  38- 
The  inflow  ratio  in  the  slipstream  is  determined  as  follows: 


Thus, 


/Ape 

=  'A*  >  -V-  C  '  —  ^ 


A-1 


^  ^  A-a 

This  expression  is  plotted  in  Figure  30* 

The  derivatives  were  given  in  Reference  12.  The  expressions 

for  these  derivatives  will  'be  summarized  here  suad  are  plotted  in 
Figures  33-37. 

Using  equation  (5)  and  noting  that  (See  Figure  l), 

-V-  A-3 


(where  n/r  here  is  the  resultant  velocity  at  the  rotor  plane  sind 
not  at  the  wing  center  of  pressure),  we  easily  obtain  the  following 
expression  for  ^  : 


a-4 
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It  is  desirable  to 


Since  this  equation  is  of  fourth  degree  in  ^ 
plot  it  in  such  a  manner  that  ^  can  be  determined  if 
and  are  known.  Dividing  the  equation  by  >0“ 


gives, 


where 


\J  CCO^O< 

>*  >  7a^=  >‘V>o  ^  >o=-f^ 


This  expression  is  plotted  in  Figure  37* 

From  the  definition  of  we  obtain, 

f^O(  ‘ 

This  is  plotted  in  Figure  33* 

Using  the  definition  of  ^  the  following  expressions  ceui  easily 
be  obtained  for  the  ^  derivatives: 


Ctrtac.<3 

X  i"  C 


bX  ^  JUA^O*  C  X^^  +  M^  cco^c<^  -V- 
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Using  the  charts  for  these  derivatives  the  induced  velocity- 
derivatives  can  easily  be  formed  from  the  following  expressions: 


1 

^4^ 


^  AxAn  CX  — 


do< 


)I 


These  are  plotted  in  Figures  3I  and-  32.  The  derivative 

from  the  chart  for  (Figure  3^). 

oGp 


differs  from 


by  only  a  constant  and  can  thus  be  found 


B.  Values  of  physical  parameters  for  Vertol  'jS  (Full  scale;  model 
scale  factor:  5.2) 


N 

=  2 

ilp 

=  . 572  ft . 

=  .075  ft. 

JU 

=  2.593  ft. 

=  1.17  ft. 

P 

=  0.3025  ft. 

<2. 

=  0.238  ft. 

Ms 

=  1.08  slug 

P 

=  .00237^  ^ 

SLR 

=  702  ft/sec 

Cw 

=  4.75  ft. 

Oo 

=  5.73 

(T 

=  0.218 

>S  =  90.^  ft  (wing  area  in  slipstream) 


A-10 

A-11 
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4 

Sw 

1 

2 

=  24.3  ft  (wing  area  outside  slipstream) 

UARe 

=  9.88 

Wt 

=  4.29  X  10^ 

=  1 

Q^qC 

=  -  0.1 

CXw 

II 

\s\ 

Kmcf 

=  1.6  X  10^ 

R 

=  4.75  ft. 

,Sr 

=  20  ft^  (area  of  horizontal  tail) 

Sf 

2 

=  20  ft  (fuselage  equivalent  flat  plate  area) 

w 

=  3160  lbs. 

S 

=  0.0087 

1 

=  12.5  ft. 

W 

=  5.^  ft. 

a 

Op 

=  3  ft.  (tail  chord) 

‘^T 

=  2° 

=  1872  (average  over  all  i-w  veQ.ues) 

•X^  nnaL)«. 

=  1885 

iTitO 

=  1858 

(main  rotor)  =  3*28 

! 

(tail  rotor)  =  2.0 

mass  of  tail  rotor  blade  =  .005  slug 

-TL 

(tail  rotor)  =  4o6  rad/sec 

Rt 

=  1.0  ft.  (radius  of  tail  rotor) 

Cpr 

=  .25  ft.  (tail  rotor  chord) 

If 

mass  main  rotor  blade  =  0.455  slug 

i4o 

-  3  •!  (tail  lift  curve  slope) 

^RT  =  5  (tail  rotor  blade  lift  curve  slope) 

C •  Determination  of  Value  for  K  from  Experimental  Data 

A  tilt-ving  VTOL  aircraft  will  not  hover  with  the  rotor  exactly 
vertical  since  the  lift  force  on  the  wing  due  to  the  slipstream 
produces  a  force  acting  in  the  negative  X-direction.  Thus  for  hover¬ 
ing  the  wing  must  be  tilted  forward  by  an  amount  such  that  the  com¬ 
ponent  of  the  wing  lift  in  the  X-directlon.  Since  the  wing  lift  is  a 
function  of  the  arbitrary  peirameter  K,  K  can  be  determined  if  it  is 
loiown  at  what  wing  tilt  angle  the  full  scale  aircraft  or  a  model  of  it 
will  hover. 

When  the  aircraft  is  hovering,  the  wing  drag  is  parallel  to  the 
thrust  and  the  lift  is  perpendicular  to  it.  Thus  we  have, 

T  <p(o  =  U  cjb^o/o  D  o(o  c-i 


where  o  —  A-<w  o  is  the  hovering  rotor  angle  of 

attack. 

Expressing  T,  L,  and  D  in  terms  of  coefficients 

p  N  TTR^  C  Ot  /IMH  P<  o 

*  -^pS  Cus  -Jf-p  5  Cps^LMndo 

as  is  given  by, 

Cus  C!uso  C 


C-2 


l4l 


C-3 


Substituting  this  into  the  above  expression  and  solving  for  K  we 


obtain 

\- 


'tojn  C  90— 


C-4 


Using  values  from  appendix  B  and  D  gives, 

3.57  ■fcaJnC^O-^^ 

The  model  of  the  Vertol  76  tested  at  the  Princeton  University 

Forward  Flight  Facility  was  found  to  hover  at  or 

o 

Cp(o~  ^  •  Thus, 

+<=  3.S7  (  .osa)  =  .'g\4 


lh2 


O  .5  1.0  IS  51.0  ^’5 
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Figure  32 
1^5 


Figure  38 
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Appendix  D 
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